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ABSTRACT 


This  thesis  supports  computer-aided  data  analysis  of  holograms  produced  from 
rocket  motor  firings.  The  work  reported  in  this  thesis  modified  existing  software  code  to 
make  it  compatible  with  installed  upgrades  in  the  microcomputer  imaging  system.  In 
particular,  this  involved  converting  the  format  of  C  language  function  calls  to  ITEX/PC 
image  processing  software  to  that  dictated  by  ITEX/PCp/us  software.  Additional  modi¬ 
fications  were  performed  to  enhance  code  portability  and  optimization.  Results  indicate 
that  all  efforts  to  incorporate  the  new  system  software  have  been  successful. 
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I.  INTRODUCTION 


A.  BACKGROUND 

This  thesis  was  a  continuation  of  the  work  currently  being  performed  at  the 
Naval  Postgraduate  School  and  was  supported  by  the  contributions  attained  in  Ref¬ 
erences  1,  2,  3  and  4.  The  goal  of  this  research  was  to  obtain  data  concerning  the 
effect  of  propellant  properties  upon  the  performance  of  solid  propellant  rocket  mo¬ 
tors.  Data  analysis  is  performed  on  holographic  images  obtained  from  combustion 
products  created  during  rocket  firings.  Holograms  are  obtained  through  the  use  of 
a  pulsed  ruby  laser  and  a  holocamera.  Reconstruction  of  the  holographic  image  for 
analysis  is  achieved  with  a  krypton  laser,  variable- power  microscope,  and  video  cam¬ 
era.  1  nfortunately.  the  optical  difTuser  required  to  record  the  hologram  creates  a 
speckle- corrupted  image  which  must  be  filtered  before  the  data  retrieval  process  is 
initiated.  The  data  attained  from  these  procedures  is  used  in  performance  studies 
stability  analysis  predictions,  and  exhaust  signature  diagnosis. 

B.  THESIS  OBJECTIVES 

The  thesis  objectives  were: 

1.  To  implement  software  modifications  in  existing  local  code  in  support  of  the 
new  computer  system  configuration. 

2.  To  enhance  program  code  portability. 

3.  To  optimize  program  code  and  attempt  to  reduce  user  interaction  where  possi¬ 
ble. 

4.  To  create  a  menu  format  for  executable  files. 

All  objectives  were  accomplished.  1TEX/PC plus  library  routines  have  been  suc¬ 
cessfully  integrated  into  the  image  analysis  programs.  The  program  file  gcnfunc.c  was 


used  to  optimize  program  size,  and  new  functions  were  created  for  user  convenience. 
A  menu  format  was  devised  to  make  the  image  analysis  procedure  easier  to  use. 


C.  SOFTWARE  SUPPORT 

The  associated  system  software  consists  of  the  following  items: 

•  Image  analysis  programs  Treated  via  thesis  research  in  References  1,  2,  3,  and 
4, 

•  Imaging  Technology’s  ImageActionplus  Software  package, 

•  Imaging  Technology’s  ITEX/PCp/tts  software  library, 

•  IBM  Disk  Operating  System  version  3.3  (DOS), 

•  Microsoft  C  optimizing  compiler  version  5.0  (MSC),  and 

•  Microsoft  Codeview  version  2.10. 

The  Image Actionp/us  software  package  is  a  menu-driven  system  that  allows  the 
user  to  process  stored  or  live  video  images  [Ref.  5].  The  two-level  menu  organization 
provides  image  analysis,  graphics,  geometric  operations,  and  image  processing  capa¬ 
bilities.  This  software  package  also  has  the  ability  to  utilize  user-generated  special 
purpose  script  files.  The  script  files  are  created  in  C  language  and  are  supported  by 
libraries  created  by  Imaging  Technology,  Inc.  Imaging  Technology’s  ITEX/PCp/us 
software  libraries  are  compatible  only  with  the  Microsoft  PASCAL,  FORTRAN,  and 
C  compilers.  These  libraries  consist  of  image  processing  and  graphics  function  sub¬ 
routines  used  in  conjunction  with  the  PCVISION plus  Frame  Grabber  board,  a  frame 
memory  and  video  digitizer  installed  in  the  IBM  PC/AT. 

D.  THESIS  ORGANIZATION 

This  thesis  is  organized  into  five  chapters.  The  introductory  chapter  provides 
the  project  background  and  goals.  It  also  discusses  thesis  objectives  and  organization. 
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Chapter  II  describes  the  system  hardware  modifications.  The  PCVISION  and 
the  PCVISION  plus  hardware  configurations  are  compared,  and  the  respective  hard¬ 
ware  capabilities  are  presented.  Also  the  hardware  performance  is  briefly  discussed. 

The  third  chapter  gives  a  synopsis  of  the  image  processing  procedure  and  dis¬ 
cusses  actual  software  programs  and  improvements  at  length.  The  chapter  also  de¬ 
scribes  the  techniques  and  software  development  tools  utilized  for  program  modifica¬ 
tion.  Finally,  Chapter  III  provides  program  specifics. 

The  software  performance  analysis  is  provided  in  Chapter  IV.  Program  analysis 
consisted  of  comparisons  in  the  following  areas:  executable  program  size,  program 
execution  times,  and  program  result  verification. 

Lastly,  Chapter  V  summarizes  the  findings  obtained  during  the  course  of  this 
research  and  provides  recommendations  for  future  thesis  research. 


II.  SYSTEM  HARDWARE  MODIFICATIONS 

A.  HARDWARE  CONFIGURATION 

Digital  image  processing  involves  the  achievement  of  the  following  objectives 
through  a  computer  interface:  recognition,  segmentation,  enhancement,  or  analysis 
of  an  image.  The  image  processing  and  digitization  system  consists  of  the  following 
components: 

•  An  IBM  PC/AT  computer  with  40  megabyte  hard  drive, 

•  An  Intel  INBOARD  386/AT  microprocessor  and  80387  math  coprocessor, 

•  An  AST  ADVANTAGE  memory  expansion  board, 

•  A  PCVISlONp/us  Frame  Grabber  board, 

•  A  computer  monitor  with  EGA  display, 

•  A  Panasonic  color  video  data  monitor,  and 

•  A  video  cassette  recorder. 

A  block  diagram  of  the  basic  elements  of  the  image  processing  system  is  presented 
in  Figure  2.1.  One  of  the  primary  components  of  this  image  processing  system  is 
the  video  cassette  recorder,  which  provides  the  introduction  of  the  image  into  the 
PCVISIONp/us  Frame  Grabber  board.  This  image  signal  consists  of  timing  and 
analog  video  information  in  a  standard  RSI 70  video  format.  Through  the  digitization 
process,  the  analog  signal  is  converted  to  a  digital  format  for  processing  and  storage. 
The  IBM  PC/AT  central  processing  unit  accesses  these  digital  values  via  the  frame 
memory. 

The  PCVISIONp/us  Frame  Grabber  is  a  single  board  that  is  connected  directly 
into  one  of  the  16-bit  expansion  slots  of  the  IBM  PC/AT  computer.  Pixels  are  stored 
in  frame  memory  after  the  Frame  Grabber  digitizes  the  analog  video  signal  at  a 
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x,y=0,0  511,0  0,0  511,0 


Figure  2.2:  Display  of  AVW  and  Hidden  Areas  of  Dual  Memory  Configu¬ 
ration.  (From  Ref.  6) 

rate  of  thirty  frames  per  second  with  eight  bits  of  accuracy.  In  previous  work,  a 
PCVISION  Frame  Grabber  board  was  utilized  (versus  the  PCVISIONp/us  version). 
The  PCVISIONp/us  board  permits  either  two  512  x  512  pixel  images  or  one  640  x 
512  pixel  image  to  be  stored  within  frame  memory  (Ref.  7:p.  1-4].  In  the  PCVISION 
board,  only  one  512  x  512  pixel  image  can  be  stored  (Ref.  8:p.  1-2].  The  bottom 
thirty-two  lines  of  the  video  output  from  both  versions  of  the  Frame  Grabber  board 
are  within  a  “hidden  area”  which  can  be  accessed  by  the  host  computer  but  not  viewed 
on  the  video  monitor  without  scrolling  the  image.  The  active  video  window  (AVW) 
is  the  area  within  frame  memory  in  which  the  pixel  values  are  displayed  or  stored  via 
the  scanning  circuit  (Ref.  6:  p.  1-14],  Figure  2.2  and  Figure  2.3  indicate  respectively 
the  dual  and  single  memory  configuration  hidden  and  AVW  areas.  Currently,  only 
the  dual  memory  configuration  is  utilized  in  thesis  research. 
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Figure  2.3:  Display  of  AVW  and  Hidden  Areas  of  Single  Memory  Config¬ 
uration.  (From  Ref.  6) 

Prior  to  display  of  the  image  on  the  video  monitor,  look-up  tables  (LUT)  are 
employed  for  the  transformation  of  pixel  intensities.  The  original  PCVISION  Frame 
Grabber  had  only  four  LUTs  associated  with  output  [Ref.  8:p.  2-8].  The  current 
Frame  Grabber  board  has  one  input  digitization  path  and  three  output  channels,  each 
allocated  eight  LUTs  [Ref.  7:p.  1-4].  These  LUTs  allow  for  real  time  processing  by 
virtue  of  their  ability  to  afford  simple  point  transformations  of  the  256  gray  levels 
without  any  processing  delay.  Look  up  tables  are  addressed  by  eight  bits  of  pixel 
information  provided  by  the  frame  memory.  A  digital-to-analog  converter  and  the 
look  up  tables  constitute  the  display  logic  unit  of  both  versions  of  the  Frame  Grabber 
boards.  The  digital-to-analog  converter  receives  gray  intensity  values  corresponding 
to  the  actual  value  addressed  within  the  LUTs  [Ref.  8:p.  2-9]. 


Both  versions  of  the  Frame  Grabber  board  are  compatible  with  the  IBM  PC/AT 
computer  bus  structure.  This  host  interface  provides  the  communication  path  be¬ 
tween  the  Frame  Grabber  boards  and  the  personal  computer.  The  previous  PC- 
VISION  board  host  interface  consisted  of  three  components:  the  Interrupt  Logic, 
Control  Register  Interface,  and  the  Frame  Memory  Interface  [Ref.  8:p.  2-1].  The 
PCVISIONp/us  host  interface  is  composed  of  the  Control  Register  Interface  and  the 
Frame  Memory  Interface  [Ref.  7:p.  4-2],  The  Interrupt  Logic  had  afforded  the 
PCVISION  Frame  Grabber  the  ability  to  be  controlled  in  real  time  while  the  IBM 
personal  computer  processed  other  instructions  [Ref.  8:p.  2-3].  This  outdated  design 
feature  has  been  superceded  by  transparent  computer  access  to  frame  memory  in  the 
PCVISION  plus  system  [Ref.  7:p.  1-6]. 

Within  the  input/output  (I/O)  space  of  the  PC/AT,  sixteen  bytes  are  reserved 
for  the  PCVISIONp/us  Frame  Grabber  control  registers.  Currently,  only  twelve  of 
these  sixteen  bytes  are  being  used  by  the  system.  Access  to  the  remaining  bytes 
is  prohibited;  they  are  reserved  for  future  system  growth  [Ref.  7:p.  4-3].  On  the 
contrary,  the  PCVISION  board  requires  32  bytes  within  the  I/O  space  for  its  control 
registers.  Only  the  first  seven  of  these  bytes  are  actually  utilized  for  control  functions 
[Ref.  8:p.  2-3]. 

The  frame  memory  of  the  PCVISIONp/us  system  is  apportioned  into  eight  equiv¬ 
alent  blocks  of  64 K  bytes.  Block  dimensions  are  512  pixels  by  128  pixels,  and  the 
physical  layout  of  the  dual-store  configuration  frame  memory  is  shown  in  Figure  2.4. 
This  configuration  can  be  contrasted  to  that  of  the  640x512  image  in  which  memory 
is  apportioned  into  eight  blocks  of  1024  pixels  by  64  pixels.  The  single-store  frame 
memory  configuration  can  be  viewed  in  Figure  2.5.  Individual  selection  of  these  blocks 
is  controlled  by  select  bits  within  the  Control  Register.  As  viewed  by  the  host,  each 
block  gives  the  semblance  consecutive  data  bytes  originating  at  the  memory  base 
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Figure  2.4:  PCVISIONp/us  Dual  Frame  Memory  Layout.  (From  Ref.  7) 


address  [Ref.  7:p.4-3).  For  the  PCVISION  system,  the  frame  memory  was  subdivided 
into  four  blocks  of  256  pixels  by  256  pixels,  forming  a  quadrant  structure  as  displayed 
in  Figure  2.6.  At  any  time  only  one  block  can  be  accessed  from  the  personal  computer 
bus  [Ref.  8:p.  2-7], 

Both  versions  of  the  Frame  Grabber  board  require  the  assignment  of  a  single 
block  of  64K  bytes  within  the  IBM  PC/AT  memory  space  located  above  640K.  Spe¬ 
cialized  software  controls  the  mapping  of  individual  blocks  of  frame  memory  into  this 
reserved  space  [Refs.  7  and  8j.  In  the  case  of  the  PCVISIONp/ws  board,  the  memory 
base  address  is  factory-configured  at  address  A0000II.  This  address  may  be  modified 
upon  board  installation  provided  the  subject  address  selected  is  w’ithin  the  guidelines 
set  forth  in  Reference  7. 

It  was  originally  intended  that  this  thesis  work  would  be  implemented  on  a 
COMPAQ  DESKPRO  3S6/20  computer  system.  The  PCVISIONp/ws  memory  base 
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re  2.5:  PCVISIONp/ws  Single  Frame  Memory  Layout.  (From  Ref.  7) 


Figure  2.6:  PCVISION  Frame  Memory  Layout.  (From  Ref.  8) 


address  was  configured  to  an  approved  value  of  D0000H.  Therefore,  the  address  region 
from  D0000H  to  DFFFFH  was  reserved  within  the  computer  memory  space,  and  all 
other  memory  and  peripherals  were  restricted  access  to  that  memory  region.  Upon 
subsequent  testing  and  diagnostic  checks,  it  was  found  that  there  was  an  unresolved 
memory  conflict  within  the  COMPAQ  computer  memory  which  affected  some  of  the 
operations  of  the  Frame  Grabber  board.  Due  to  necessity,  the  board  was  consequently 
installed  in  the  IBM  PC/AT.  Memory  base  address  D0000H  was  again  selected;  this 
time  there  were  no  resident  memory  conflicts. 

B.  HARDWARE  PERFORMANCE 

With  the  introduction  of  the  PCVISIONp/us  Frame  Grabber  board,  the  ability 
to  process  two  separate  images  in  frame  memory  is  invoked.  The  dual  memory  con¬ 
figuration  allows  the  host  to  select  either  image  area  for  analysis  and  manipulation. 
A  display  memory  function  allows  the  user  to  view  independently  either  image  area. 
These  capabilities  allow  for  considerable  future  flexibility  in  the  analysis  of  the  rocket 
motor  holograms. 

The  PCVISIONp/us  board  not  only  provides  more  image  processing  capabilities 
than  the  PCVISION  board,  but  it  also  establishes  the  need  for  additional  function 
libraries  to  support  its  abilities.  For  this  reason,  image  analysis  programs  devised  in 
prior  thesis  research  had  to  be  modified  to  reflect  the  new  supporting  libraries.  These 
programs  and  their  resulting  modifications  are  emphasized  in  the  following  chapter. 
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III.  SOFTWARE  PROGRAMS  AND 
IMPROVEMENTS 


A.  IMAGE  PROCESSING  PROCEDURE 

As  was  discussed  in  the  introductory  chapter,  holographic  techniques  are  em¬ 
ployed  to  obtain  the  image  of  combustion  products  from  small  rocket  motors  during 
firing.  During  the  hologram  reconstruction  process,  a  krypton  laser  is  used  in  con¬ 
junction  with  a  variable  power  microscope.  Due  to  the  presence  of  a  diffuser  that  is 
required  to  eliminate  thermal  gradient  effects,  speckle  noise  is  introduced  into  the  cap¬ 
tured  image.  The  following  paragraphs  describe  the  procedure  developed  to  acquire 
and  process  this  image  for  statistical  study. 

1.  Image  Acquisition 

The  speckle-corrupted  image  is  recorded  for  future  processing  via  a  combi¬ 
nation  of  a  video  cassette  recorder  and  low-light-level  camera  attached  to  the  variable 
power  microscope.  The  magnification  power  of  the  microscope  determines  the  rel¬ 
ative  size  of  the  pixel  elements  comprising  the  image.  Pixels  are  the  standard  unit 
for  the  measurement  of  the  features  within  the  image.  Object  sizes  are  determined 
by  scaling  the  measured  object  width  (in  pixels)  by  the  dimension  spanned  by  one 
pixel  (in  microns  per  pixel).  A  quantization  error  is  introduced  into  this  calculation 
due  to  the  necessity  of  measuring  dimensions  with  an  integer  number  of  pixels.  As 
the  width  decreases,  the  proportional  level  of  uncertainty  increases,  which  produces 
a  more  erroneous  scale  factor  calculation.  A  calibration  object  is  inserted  within  the 
holographic  image  to  determine  the  width  spanned  by  a  pixel.  It  is  essential  that  this 
object  covers  an  appreciable  portion  of  the  monitor  screen  (at  least  one  half  of  the 


area)  to  counteract  this  quantization  error.  The  calculation  of  the  conversion  factor 
from  pixels  to  microns  becomes  increasingly  important  in  subsequent  steps  of  image 
processing. 

2.  Image  Digitization 

From  the  VCR  tape  obtained  from  the  hologram  reconstruction,  the  menu- 
driven  ImageActionp/us  software  is  used  to  “grab”  and  digitize  the  desired  image. 
The  PCVISIONp/us  Frame  Grabber  board  implements  this  procedure.  The  digitized 
video  image  is  comprised  of  a  512  x  480  array  of  pixels  having  gray  intensity  levels 
ranging  from  0  to  255.  These  intensity  levels  extend  from  blackest-black,  level  0,  to 
whitest-white,  level  255. 

3.  Speckle  Reduction  Filtering 

Image  filtering  is  the  paramount  technique  used  in  speckle  reduction.  The 
amount  of  speckle  reduction  achieved  by  filtering  algorithms  is  gauged  by  a  calculation 
of  the  speckle  index.  Speckle  index  is  defined  as  the  ratio  of  the  image’s  local  deviation 
to  the  local  mean  of  pixel  values,  and  this  measure  is  indicative  of  the  random  speckle 
noise  inherent  in  the  image  [Ref.  l:pp.  17—19). 

The  three  filtering  algorithms  are  based  upon  synthetic  aperture  radar 
(SAR)  theory  presented  in  References  9  and  10.  There  are  three  separate  filter  mod¬ 
els.  Derived  from  nonlinear  techniques,  two  filter  designs  are  based  upon  statistical 
methods.  The  third  filter  model  has  its  foundation  in  a  geometric  hulling  algorithm. 
The  geometric  hulling  algorithm  produces  the  best  results  by  enabling  the  user  to 
discontinue  filtering  at  a  less  severe  level  of  resolution  degradation. 

4.  Image  Threshold 

In  an  ideal  image,  feature  particles  and  the  background  would  have  distinct 
gray  intensity  levels.  The  process  of  thresholding  creates  a  binary  image  in  which 
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image  features  will  appear  black  and  the  background  white.  In  reality,  the  thresholded 
image  will  contain  some  amount  of  speckle  manifested  as  feature  particles  and  will 
have  lost  portions  of  actual  particles  during  the  filtering  and  thresholding  process. 
Image  evaluation  errors  are  inherent  whenever  the  grey  levels  of  the  object  overlap 
the  gray  levels  of  the  image  background.  The  selection  of  a  threshold  value  is  obtained 
from  observation  of  the  filtered  image.  ITEX/PCp/us  software  allows  the  user  to  select 
a  threshold  value  for  which  the  designated  LUT  is  modified.  This  modification  is  not 
made  permanent  until  the  image  is  mapped  into  frame  memory,  allowing  the  user 
to  rapidly  iterate  values  to  make  a  determination  on  the  most  appropriate  threshold 
value. 

5.  Feature  Identification 

Once  the  filtered  image  has  been  thresholded,  it  is  ready  to  undergo  fea¬ 
ture  identification.  It  is  imperative  that  the  programmed  image  has  been  previously 
thresholded  due  to  the  logic  of  the  identification  algorithm.  This  procedure  locates 
all  connected  adjacent  feature  pixels  and  assigns  feature  value  numbers  to  them. 
Through  this  labeling  process,  interconnected  pixels  are  identified  as  an  individual 
object.  The  total  number  of  objects  recognized  serves  as  input  for  the  next  step 
of  image  processing.  This  number  is  restricted  to  32,767  or  below  based  upon  the 
limitations  imposed  by  the  C  language  on  integer  values  [Ref.  4:p.  9]. 

6.  Feature  Sizing 

The  procedure  of  feature  sizing  determines  both  the  maximum  vert  ical  and 
horizontal  diameters  of  the  identified  objects.  The  total  area  of  the  subject  particles 
is  computed,  and  all  these  measurements  are  delineated  in  terms  of  microns.  As 
described  before,  a  conversion  factor  is  necessary  for  manipulation  of  the  feature 
particle's  pixel  measurements  into  actual  physical  dimensions.  This  constant,  called 
SCALE-FACTOR,  is  presently  based  upon  the  dimensions  of  a  threaded  screw  at  the 
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same  magnification  of  the  captured  image.  The  screw  threads  serve  as  a  calibration 
object  to  transform  image  object  pixel  size  data  into  actual  physical  measurement 
criteria.  The  final  output  of  the  feature  sizing  routine  is  presented  in  tabular  form 
suitable  for  input  into  the  STATGRAPHICS  software  package,  a  statistical  analysis 
tool. 

7.  Histogram  production 

To  exhibit  feature  size  data  in  a  histogram  form  for  ready  evaluation,  the 
STATGRAPHICS  software  package  is  utilized.  This  particular  software  package  was 
selected  because  it  is  compatible  with  the  IBM-PC/AT  system  and  is  relatively  easy 
to  use.  Histograms  produced  by  STATGRAPHICS  are  used  for  further  statistical 
analysis  to  determine  total  particle  distributions  of  holograms. 

B.  PROGRAM  OPTIMIZATION 

The  primary  purpose  of  this  thesis  research  was  to  modify  existing  C  language 
routines  to  make  them  compatible  with  the  ITEX/PCp/us  software.  Subsequent  pro¬ 
gram  optimization  was  performed  in  support  of  this  goal.  One  significant  improve¬ 
ment  was  the  removal  of  all  general-use  functions  from  the  main  programs.  Although 
a  general  support  functions  program  file  had  been  created  in  the  previous  rendition, 
this  file  had  not  been  adequately  utilized.  As  a  result,  all  the  main  programs  had 
needlessly  large  file  sizes. 

The  C  language  provides  the  ability  to  use  the  infinitely-abusable  GO  TO  state¬ 
ment.  This  creates  a  situation  in  which  branching  occurs  in  a  deeply  nested  structure. 
The  GO  TO  organization  may  prove  effective  in  some  instances,  but  it  always  pro¬ 
duces  confusion  for  the  individual  having  the  misfortune  to  try  to  modify  that  code 
sequence.  For  this  reason,  program  structures  were  modified  to  eliminate  all  GO  TO 
branch  statements,  thereby  making  code  more  modular  and  portable. 
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A  menu  format  was  developed  for  use  with  the  executable  programs.  This 
format  will  facilitate  demonstrations  of  image  processing  and  makes  the  routines 
“user- friendly” .  In  addition,  supporting  functions  were  devised  for  some  of  the  pro¬ 
grams  to  further  eliminate  the  need  for  user  interaction  with  the  programs  once  the 
selected  image  had  been  loaded  and  image  processing  regions  were  identified.  These 
modifications  facilitated  more  rapid  program  execution.  A  more  detailed  discussion 
of  program  modification  specifics  will  occur  in  a  succeeding  section. 

C.  SOFTWARE  DEVELOPMENT  TOOLS 

The  Microsoft  Codeview  Debugger  program  was  utilized  for  program  develop¬ 
ment.  It  requires  that  specific  compiler  directives  be  used  during  program  compi¬ 
lation.  Although  Codeview  is  an  excellent  developmental  tool,  program  debugging 
can  become  a  tedious  process.  No  program  modifications  can  be  performed  while  in 
Codeview,  and  for  that  reason,  it  is  not  an  adequate  tool  for  what  is  needed  for  pro¬ 
gram  modules  of  this  size.  The  ITEX/PC plus  routines  require  large  memory  module 
handlers,  and  few  other  debuggers  are  able  to  meet  this  requirement.  Quick  C  version 
2.0  is  capable  of  handling  large  memory  modules  and  has  an  online  source  debugger 
which  is  more  versatile  than  Codeview  [Ref.  11).  It  is  recommended  that  in  future 
work  with  these  program  modules,  Quick  C  should  be  included  in  system  upgrades. 

Several  batch  files  were  created  for  compilation  and  linker  processes  since  Code- 
view  required  special  compiler  directives.  There  are  two  different  versions  of  these 
files:  compilation  and  link  with  debugger  information,  and  compilation  without  de¬ 
bugger  information  for  file  optimization  and  compaction.  Listings  and  justifications 
for  these  DOS  Batch  files  are  referenced  in  Appendix  K. 
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D.  PROGRAM  SPECIFICS 

The  succeeding  paragraphs  discuss  the  programs  that  were  developed  in  support 
of  this  thesis  research.  These  existing  programs  were  modified  to  be  compatible  with 
the  ITEX/PCp/us  software  libraries  and  the  Image  Actionp/us  routines.  Where  at 
all  possible,  program  optimization  was  performed.  In  some  cases,  code  had  to  be 
completely  revamped  to  accommodate  new  requirements.  All  source  code  listings  are 
included  as  appendices  at  the  end  of  this  document.  A  block  flow  diagram  of  program 
usage  is  provided  in  Figure  3.1. 

1.  Program  File  thesis. h 

The  program  header  file  thesis. h  (Appendix  A)  contains  all  requisite  library 
include  files,  C  manifest  constants,  and  MACRO  definitions  as  necessitated  by  main 
programs.  For  code  portability,  it  is  imperative  that  a  centralized  copy  of  all  common 
declarations  and  definitions  be  maintained  as  the  programs  evolve.  Therefore,  any 
modifications  to  common  parameters  can  be  performed  by  substitutions  within  this 
file. 

With  the  advent  of  ITEX/PCp/us  software,  additional  include  files  were 
necessary  beyond  the  scope  of  the  original  header  file  listed  in  Reference  4.  For 
example,  ITEX/PC  libraries  require  the  use  of  the  include  file  itexpc.h  [Ref.  12:p. 
1-15].  This  requirement  has  been  superceded  by  the  ITEX/PCp/us  software  which 
dictates  the  use  of  itexpfg.h  and  stdtyp.h  include  files.  Not  only  must  these  include 
files  be  listed  in  the  header  file,  but  to  be  utilized,  they  must  also  be  resident  in  the 
actual  program  directory  [Ref.  6:p.  1-9]. 

The  added  capability  of  dual  memory  required  that  many  of  the  calls  to  the 
software  libraries  be  altered.  Additional  manifest  constants  were  needed  for  memory 
specification.  They  included  MEMORY,  defining  the  type  of  frame  memory  con¬ 
figuration.  and  constants  necessary  to  define  the  size  of  the  selected  frame  memory 
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configuration:  XSIZE,  YSIZE,  and  DEPTH.  XSIZE  and  YSIZE  indicate  respectively 
the  horizontal  and  vertical  coordinates  of  the  memory  size,  and  DEPTH  refers  to  the 
number  of  bits  assigned  per  pixel  in  frame  memory  [Ref.  6:p.  2-6]. 

2.  Program  File  geofil.c 

Appendix  B  contains  the  listing  for  geofil.c ,  a  program  based  upon  a  geo¬ 
metric  hulling  algorithm  introduced  in  References  9  and  10.  The  algorithm  is  applied 
to  the  image  and  its  complement  by  performing  comparisons  between  a  central  pixel 
value  and  adjacent  pixels  in  horizontal,  vertical  and  diagonal  directions.  The  original 
geometric  filter  routines  were  adapted  by  Edwards  [Ref.  1]  from  information  set  forth 
in  References  9  and  10;  and  Kaeser  successfully  adapted  these  to  the  C  language  [Ref. 
4].  Edwards’  FORTRAN  filtering  routine  was  capable  of  handling  only  a  quarter  of 
the  actual  image  due  to  memory  restrictions,  whereas  the  C  language  version  fully 
supported  the  entire  image.  A  large  part  of  the  C  version’s  success  was  due  to  the 
utilization  of  virtual  memory  arrays  for  data  processing.  These  arrays  are  stored  on 
disk,  but  they  are  accessed  as  though  they  are  stored  in  memory.  File  constraints 
within  the  associated  operating  system  pose  the  only  limitation  upon  the  virtual  array 
size  [Ref.  13]. 

Few  modifications  were  performed  on  the  geofil.c  program  to  make  it  com¬ 
patible  with  the  ITEX/PC plus  software.  The  majority  of  the  changes  for  compatibil¬ 
ity  encompassed  general  support  functions.  However,  much  of  the  C  language  code 
was  altered  to  attain  optimum  modularity.  This  filter  program  had  incorporated  the 
use  of  GO  TO  statements  to  facilitate  branching  between  the  filter  stepping  functions. 
Alterations  included  the  removal  of  these  subject  statements  and  the  formation  of  two 
recursive  functions.  These  recursive  functions  support  program  modularity  and  per¬ 
mit  easier  modifications  for  future  updates. 
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One  common  modification  performed  on  geofil.c  and  its  other  main  pro¬ 
gram  counterparts  involved  the  declaration  of  the  image  filename  and  comment  line 
(comline).  The  header  file  used  by  ITEX/PC,  itexpc.fi,  required  special  file  type  struc¬ 
tures  used  in  conjunction  with  parameter  passing  for  SAVIM  and  READIM  functions. 
These  structures  consisted  of  strings  with  length  20  or  200  characters,  and  were  clas¬ 
sified  as  typedef  LS20  or  LS200,  respectively  [Ref.  12:p.  1-19).  ITEX/PCp/us  limited 
filename  to  20  characters  and  comment  lines  to  200  characters,  but  removed  the 
special  structure  requirement.  SAVIM  and  READIM  routines  now  require  that  the 
image  filename  and  comment  line  be  classified  as  pointers  of  type  character  [Ref.  6]. 
Consequently,  all  program  files  were  altered  to  reflect  this  requirement. 

3.  Program  File  Istat.c 

The  program  Istat.c  (Appendix  C)  uses  local  statistical  methods  to  calcu¬ 
late  a  weight,  k,  which  determines  the  gray  level  of  a  pixel  [Refs.  9  and  10].  The 
mean  and  variance  local  estimates  within  a  5x5  window  encompassing  the  selected 
pixel  determine  the  value  of  k.  The  k  value  then  specifies  the  new  gray  level  of  the 
pixel  in  question  [Ref.  l:p.  24],  Virtual  arrays  again  assist  in  the  data  manipulation 
and  storage  processes. 

Kaeser’s  version  of  this  filtering  algorithm  required  that  the  user  input  the 
standard  deviation  of  the  image  directly  into  the  program  [Ref.  4].  The  operator  had 
to  obtain  this  value  from  the  histogram  function  of  the  ImageAction  software  [Ref. 
14:p.  11-9-3]  before  running  the  filter  program.  ImageActionp/us  software  no  longer 
provides  the  image  mean  and  standard  deviation  calculations  [Ref.  5].  Consequently, 
a  new  function,  dcv(stddcv).  was  devised  to  calculate  the  image  standard  deviation 
and  input  it  directly  into  the  filter  routine  without  user  intervention.  This  new 
function  presents  a  vast  improvement  in  terms  of  user  convenience. 
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4.  Program  File  2sigma.c 

The  2sigma  filter  algorithm  provides  the  basis  for  2sigma.c  (Appendix  D) 
and  was  adapted  from  the  program  introduced  by  Edwards  [Ref.  1].  From  the 
standard  deviation  of  a  Gaussian  distribution,  it  is  known  that  95.5%  of  a  selected 
region’s  pixels  will  be  within  two  standard  deviations  of  the  local  mean.  An  averaging 
formula  has  been  devised  in  which  pixels  meeting  the  “2sigma”  criteria  are  averaged 
for  a  given  pixel,  and  those  outside  the  range  are  rejected  [Ref.  10].  These  calculations 
are  performed  for  5x5  arrays  of  pixels.  The  program  2sigma.c  utilizes  virtual  arrays 
and  the  standard  deviation  function  as  well.  It  also  has  reduced  size  due  to  the 
removal  of  all  general  supporting  functions  into  the  program  file  genfunc.c. 

5.  Program  File  threshit.c 

Threshit.c  (Appendix  E)  produces  a  binary  version  of  the  filtered  image. 
The  program  operator  enters  a  threshold  value  which  sets  all  pixel  values  above  it  to 
WHITE  (level  255),  and  all  those  below  to  BLACK  (level  0).  The  user  then  views 
the  thresholded  image  and,  if  the  result  is  not  satisfactory,  is  given  the  option  to 
enter  a  new  threshold  limit.  Threshit.c  utilizes  a  function  named  thrcshii( )  which  is 
physically  located  within  the  genfunc.c  program  file.  This  organization  was  due  to 
the  fact  that  other  program  files,  such  as  feat. id. c ,  call  the  thrcshit()  function. 

The  key  to  the  threshold  process  is  a  function  call  in  the  ITEX/PCp/us soft¬ 
ware  package  called  THRESHOLD.  The  ITEX/PC  function  is  also  named  THRESH¬ 
OLD  and  has  the  following  format: 

void  threshold(/oi/Tuf.  highevt)-, 

The  lowcut  and  higheui  integer  values  defined  the  lower  and  upper  boundaries  of  the 
thresholded  region.  This  threshold  procedure  could  only  be  used  on  the  currently 


selected  LUT  [Ref.  12:p.  7-27].  The  function  in  ITEX/PCp/tis  software  used  to  set  a 
selected  LUT  to  a  prescribed  threshold  is: 

void  threshold(proup,  bank,  highcut,  lowcut)\ 

where, 

group  specifies  the  LUT:  RED,  GREEN,  BLUE,  or  INPUT, 
bank  specifies  the  output  LUT  used,  in  the  range  of  0  to  7, 
highcut  specifies  the  upper  bound  for  the  threshold,  and 
lowcut  specifies  the  lower  bound  for  the  threshold. 

The  threshit.c  program  utilizes  all  four  LUT  groups  and  the  Oth  256-byte  bank  to 
specify  a  linear  transform  output.  This  specification  differs  significantly  from  the 
original  ITEX/PC  function  in  the  amount  of  versatility  granted  the  programmer. 
Any  of  the  LUTs  can  be  modified  with  this  threshold  command,  not  just  the  current 
LUT  [Ref.  6:p.  7-35].  The  lowcut  value  must  never  be  specified  as  higher  than  highcut 
value,  but  these  values  may  be  made  equivalent  to  display  one  intensity  at  level  255 
(WHITE). 

Once  the  user  is  satisfied  with  a  particular  image  threshold,  the  image  can 
be  mapped  to  memory.  This  procedure  involves  the  use  of  another  ITEX/PCp/us 
function: 

int  maplut(<7ro up,  bank,  x,  y,  dx,  dy)\ 

MAPLUT  modifies  actual  memory  values  rather  than  just  LUT  values  as  THRESH¬ 
OLD  does.  The  group  and  bank  variables  are  identical  to  those  used  in  the  THRESH¬ 
OLD  function,  but  now  only  the  INPUT  LUT  is  modified  for  mapping  purposes  in 


threshit.c.  The  quantities  x,  y,  dx,  and  dy  refer  to  the  horizontal  and  vertical  coor¬ 
dinates  of  the  512x512  image  [Ref.  6:p.  7-17],  As  in  the  case  of  the  THRESHOLD 
function,  ITEX/PC  permits  the  user  to  transform  only  that  area  specified  by  the 
current  LUT.  The  variables  group  and  bank  are  not  defined  in  ITEX/PC  MAPLUT 
function  [Ref.  12:p.  7-14]. 

6.  Program  File  feat-id.  c 

It  is  mandatory  that  an  image  is  thresholded  prior  to  initiating  the  program 
feat-id. c  (Appendix  F).  For  this  reason,  the  user  is  given  the  option  within  feat-id. c 
to  threshold  the  image  if  thresholding  was  not  completed  previously.  Conversion  to 
a  feature-labeled  image  requires  an  iterative  process,  whereby  the  image  is  evaluated 
row  by  row.  Once  this  feature  identification  process  is  complete,  the  image  can  be 
further  processed  for  specific  feature  data. 

The  feat-id. c  program  necessitated  few  modifications  for  compatibility  with 
ITEX/PCp/us  software.  Again,  the  more  extensive  changes  were  performed  on  the 
general  support  functions.  Feat-id. c  used  the  ITEX/PCp/tis  functions  WPIXEL 
and  RPIXEL  extensively,  but  these  calls  retained  the  same  format  as  that  used  in 
ITEX/PC  software  [Refs.  6  and  12]. 

7.  Program  File  sizeit.c 

Program  sizeit.c  (Appendix  G)  creates  an  output  file  size.dat  which  con¬ 
tains  the  dimensions  of  all  image  features  in  a  format  necessary  for  input  into  a  sta¬ 
tistical  analysis  program.  Sizeit.c  relies  upon  pointer  structures  for  dynamic  memory 
allocation  of  processed  feature  data. 

A  conversion  factor.  SCALE-FACTOR,  figures  prominently  in  this  sizing 
algorithm.  Kaeser  determined  this  factor  to  be  3.7353  contingent  on  image  calibration 
[Ref.  4],  Upon  subsequent  runs  of  the  sizeit.c  program,  however,  a  SCALE-FACTOR 
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of  10.0000  was  necessary  to  achieve  compatible  results  with  Kaeser’s  program  [Ref. 
4:p.  17].  The  constant  term  SCALE-FACTOR  is  defined  in  the  header  file  thesis. h 
and  can  be  easily  modified  if  required.  Based  upon  equipment  configuration  during 
hologram  processing,  the  subject  conversion  factor  is  still  questionable  in  terms  of 
accuracy.  For  this  reason,  an  additional  function  was  devised  to  determine  a  con¬ 
version  factor  based  upon  user  input.  This  function  was  entitled  sca.le(),  and  its 
characteristics  are  discussed  in  the  section  relevant  to  genfunc.c. 

8.  Program  File  speckle. c 

To  calculate  the  value  of  speckle  index  (SI)  after  each  filter  iteration, 
speckle. c  (Appendix  H)  is  called  to  evaluate  filtering  effectiveness.  This  calculation 
is  performed  upon  specified  dimensions  of  the  filtered  image;  integer  values  for  im¬ 
age  row  and  column  are  passed  to  this  function  from  the  filter  routines.  Values  for 
speckle  index  typically  range  from  a  maximum  value  of  1.0,  a  rare  occurrence,  to 
a  hypothetical  value  of  0  [Ref.  l:p.  20].  No  program  modifications  were  necessary 
for  ITEX/P Cplus  compatibility.  Speckle  index  is  obtained  purely  from  statistical 
manipulations. 

9.  Program  File  vir.arry.c 

Vir.arry.c  (Appendix  I)  consists  of  a  group  of  functions  supporting  setup 
and  access  of  the  virtual  disk-based  array  stratagem  employed  primarily  by  the  filter 
routines.  The  virtual  array  functions  were  adapted  by  Kaeser  [Ref.  4]  from  Refer¬ 
ence  8.  These  arrays  reside  on  disk  and  are  accessed  as  though  they  are  located  in 
memory.  This  feature  overcomes  the  memory  limitation  normally  corresponding  to 
array  usage.  Through  the  use  of  MACROS  and  pointer  notation,  the  C  programming 
language  provides  the  means  of  managing  virtual  arrays  [Ref.  8].  Thesis. h  contains 
the  obligatory  MACROS  needed  for  structure  definition.  All  virtual  array  functions 


mandated  the  use  of  long  integers  for  indexing  purposes  [Ref.  13).  Since  the  dimen¬ 
sions  of  the  actual  image  are  defined  in  terms  of  integer  values,  a  conversion  interface 
problem  developed  in  all  the  filter  programs.  This  was  further  complicated  by  the 
fact  that  the  ITEX/PCp/us  WPIXEL  function  references  datum  points  in  terms  of 
integer  values.  This  function  is  essential  for  “writing”  the  filtered  pixel  values  to 
memory  from  the  virtual  array  routines  [Ref.  6:p.  5-4].  With  the  use  of  the  C 
language  “(long)”  command,  the  integer  indices  referencing  pixel  values  within  the 
image  can  be  used  simultaneously  to  reference  values  within  the  virtual  arrays  with 
no  conversion  difficulty  [Ref.  15:p.  36]. 

10.  Program  File  genfunc.c 

The  program  file  genfunc.c  (Appendix  J)  contains  image  processing  general 
support  functions.  These  functions  maintain  routine  evolutions  occurring  within  the 
main  image  processing  programs.  The  following  sections  discuss  functions  that  were 
either  created  or  modified  due  to  the  advent  of  the  ITEX/PCp/us  libraries. 

a.  Function  startitf) 

This  function  makes  provisions  for  initial  Frame  Grabber  board  setup. 
All  the  main  function  files  with  the  exception  of  speckle. c  and  vir.arry.c  access  this 
function.  Due  to  the  introduction  of  the  dual  memory  configuration,  new  ITEX/PCp/us 
calls  were  required  in  addition  to  what  Kaeser  utilized  [Ref.  4].  The  SETDIM  rou¬ 
tine  determines  the  configuration  of  the  memory,  either  dual-store  or  single-store  [R?f. 
6:p.  2-6].  SELECT.MEM  determines  which  frame  memory  will  be  accessed  while  in 
the  dual-store  mode.  Frame  memory  choices  consist  of  MEM.A  or  MEM.B  as  de¬ 
fined  in  the  itexpfg.h  header  file  [Ref.  6:p.  2-19].  DISPLAY _MEM  must  be  used  in 
conjunction  with  the  SELECT.MEM  routine.  It  specifies  which  memory  area  will  be 
displayed  and  gives  the  user  the  option  of  writing  to  one  image  area  while  displaying 


the  other  [Ref.  6:p.  2-20].  For  the  current  application,  MEM_A  is  both  the  accessed 
and  displayed  frame  memory. 

b.  Function  readit() 

The  image  is  read  from  disk  storage  and  placed  on  the  video  monitor 
via  the  rtadit()  function.  Any  comment  associated  with  the  image  when  it  was  stored 
will  be  displayed  at  that  time.  The  ITEX/PCp/us  function  READIM  is  the  key  to 
this  procedure.  The  READIM  call  format  is  relatively  unchanged  from  the  ITEX/PC 
version.  If  an  error  occurs  when  a  file  is  being  read,  the  function  will  return  an  error 
code  to  identify  why  the  operation  failed.  The  itcxpfg.h  header  file  defines  these  error 
code  values  [Ref.  6:p.  6-5].  The  standard  ITEX/PC  format  did  not  make  provisions 
for  these  values  in  its  associated  header  file,  itexpc.h\  therefore,  file  errors  are  referred 
to  numerically  [Ref.  12:p.  6-7].  Accordingly,  the  readit()  function  was  modified  to 
reflect  the  error  nomenclature. 

c.  Function  saveit() 

The  saveitf)  function  allows  the  user  to  save  the  image  to  disk  memory 
via  the  ITEX/PCp/us  function  SAVIM.  The  image  can  be  stored  using  compression; 
although  compression  may,  in  rare  cases,  cause  the  image  to  occupy  more  disk  space 
than  an  uncompressed  version.  If  compression  is  used,  the  file  will  be  automatically 
decompressed  when  read  from  the  disk.  As  in  the  case  of  the  READIM  function. 
SAVIM  returns  error  code  values  if  the  save  operation  is  unsuccessful.  Subsequent 
modifications  were  performed  to  the  saveit()  code  to  make  it  compatible  with  the 
new  documentation  [Refs  6  and  12]. 

d.  Function  dtv(stddtv) 

Both  program  files  Isiat.c  and  Ssigma.c  require  the  computation  of  the 
standard  deviation  of  the  image  for  their  respective  filtering  algorithms.  The  function 


dev(stddev)  computes  this  value  through  an  iterative  process  utilizing  the  RPIXEL 
function.  Since  the  standard  deviation  is  computed  directly  from  the  image  without 
user  input,  there  is  a  slight  increase  in  filter  pr  mn  execution  times  associated  with 
the  implementation  of  this  function. 

e.  Function  scale(factor) 

This  function  supports  the  sizeit.c  program  by  allowing  the  user  flex¬ 
ibility  in  the  determination  of  a  scale  factor  used  in  the  sizing  algorithm.  The  user 
is  given  the  option  of  using  the  SCALE-FACTOR  as  delineated  in  the  header  file 
thesis. h,  or  of  entering  dimensions  of  a  calibration  object.  If  the  user  decides  upon 
the  latter,  the  length  of  the  object  must  be  entered  in  terms  of  pixels  and  inches. 
The  function  will  then  calculate  the  new  SCALE.FACTOR  and  pass  this  value  to  the 
sizeit.c  program. 

E.  SUMMARY 

The  modifications  and  improvements  performed  on  the  image  processing  pro¬ 
grams  were  necessary  to  make  them  compatible  with  the  environment  established  by 
the  installment  of  the  PCVISIONp/us  Frame  Grabber  board  into  the  IBM  PC/AT. 
These  subject  programs  support  all  requirements  for  image  acquisition,  filtering, 
thresholding  and  particle  identification.  The  following  chapter  examines  the  per¬ 
formance  and  verification  of  these  programs  at  length. 


IV.  SOFTWARE  PERFORMANCE  ANALYSIS 

This  chapter  will  discuss  the  pertinent  details  concerning  software  performance. 
The  subjects  of  executable  program  size,  timing  analysis,  and  program  verification 
will  be  examined  at  length  in  the  following  sections.  Any  departures  from  results 
obtained  from  the  original  C  language  programs  will  be  appraised. 

A.  EXECUTABLE  PROGRAM  SIZE 

Table  4.1  shows  the  comparison  between  executable  program  size  of  ITEX/PC 
and  ITEX/PCp/us  supported  software.  The  size  information  is  presented  in  terms  of 
kilobytes  (KB)  and  was  acquired  via  the  use  of  the  DOS  command,  DIR  [Ref.  4]. 
The  ITEX/PC-based  programs  written  by  Kaeser  [Ref.  4]  have  significantly  smaller 
executable  file  sizes.  This  is  attributed  to  the  fact  that  during  the  LINK  process, 
Kaeser  linked  main  programs  to  required  library  files  only.  It  was  not  necessary  to 
link  to  the  object  modules  of  the  gcnfunc.c ,  speckle. e,  or  the  v ir.arry.c  files  since  all 
necessary  functions  were  appended  to  the  end  of  the  main  program  files  [Ref.  4].  In 
addition,  the  ITEX/PC  function  library  differs  from  the  ITEX/PC  plus  version,  and 
these  differences  will  also  impact  executable  file  size. 

TABLE  4.1:  Comparisons  of  Executable  Program  Size  (KB) 


FILENAME 

ITEX/PC 

ITEX/PC  plus 

THRESHIT.EXE 

31,467 

63,101 

FEATJD.EXE 

32.509 

64.407 

SIZEIT.EXE 

46,281 

65,429 

2SIGMA.EXE 

49,483 

65.491 

LSTAT.EXE 

49,475 

65.499 

GEOFIL.EXE 

49.969 

66.317 
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TABLE  4.2:  Comparisons  of  Actual  Program  Size  (KB) 


FILENAME 

ITEX/PC 

ITEX/PC  plus 

THRESHIT.C 

5,807 

1,336 

FEATJD.C 

10,130 

6,003 

SIZEIT.C 

11,011 

7,426 

2SIGMA.C 

14,623 

5,396 

LSTAT.C 

14,434 

5,396 

GEOFIL.C 

14,888 

7,916 

One  stated  objective  of  this  thesis  was  to  optimize  the  existing  C  language 
code.  Although  executable  program  size  comparisons  contradict  this  goal,  it  is  shown 
in  Table  4.2  that  actual  program  sizes  for  the  ITEX/PCp/us  compatible  programs  are 
substantially  reduced  from  the  original  versions.  Two  linker  options,  /F  and  /PAC, 
were  initially  employed  to  achieve  more  compact  executable  files.  The  use  of  these 
options  should  produce  more  rapid  execution  of  files  and  shorter  program  load  times 
[Ref.  16:p.  185].  The  / E  option  had  been  utilized  by  Kaeser  in  Reference  4,  but 
this  file  packing  option  may  not  always  give  satisfactory  results.  It  can  occasionally 
increase  file  size,  and  this  obviously  reduces  savings  in  disk  space  [Ref.  17:p.  268]. 

Upon  linking  program  files  with  the  /PAC  and  /F  options,  it  was  found  that 
these  options  did  not  pose  an  improvement  over  the  /E  option.  In  fact,  executable 
file  sizes  were  increased.  For  example,  the  gcofil.fH  file  was  73.132  KB  when  these 
options  were  used  in  the  LINK  process.  The  /E  option  produced  an  executable  file 
size  of  66.317  KB.  All  other  executable  files  were  affected  similarly.  Accordingly,  the 
/E  option  was  utilized  for  the  LINK  procedure.  Appendix  K  contains  the  listings  of 
the  batch  files  used  to  compile  and  link  the  program  files. 


TABLE  4.3:  Program  Execution  Times 


PROGRAM 

TIME(SEC) 

FEATURE  IDENTIFICATION 

809.82 

FEATURE  SIZING 

1780.03 

2SIGMA  FILTER 

760.55 

LOCAL  STATISTICS  FILTER 

1324.47 

GEOMETRIC  FILTER 

3403.13 

B.  PROGRAM  EXECUTION  TIMES 

The  jl7res4-img  library  image  was  used  to  test  program  execution  times.  This 
was  the  same  image  as  that  used  by  Kaeser  in  Reference  4,  depicting  an  Air  Force  Res¬ 
olution  Target  (AFRT)  obscured  by  speckle  noise.  The  unfiltered  image  jl7res4-ing 
is  depicted  in  Figure  4.1.  Testing  was  accomplished  by  the  use  of  a  timer  program, 
timer. c,  in  conjunction  with  a  DOS  input  file  to  eliminate  user  intervention.  Previous 
thesis  work  involved  time  testing  of  one  quarter  of  the  512x512  image.  MSFORT 
versions  of  the  programs  could  effectively  handle  only  one  quadrant  of  the  full  screen 
image.  Consequently,  Kaeser  conducted  execution  time  comparisons  on  only  one 
quarter  of  the  image  [Ref.  4].  Due  to  the  banded  nature  of  the  ITEX/PC plus  frame 
memory,  quadrant  timing  tests  of  the  image  would  not  produce  results  that  could  be 
compared  with  the  original  findings  [Ref.  5];  therefore,  timing  analysis  of  the  pro¬ 
grams  was  performed  on  the  entire  image.  The  resultant  execution  times  are  listed 
in  Table  4.3. 

C.  PROGRAM  VERIFICATION 

To  verify  the  results  obtained  from  the  ITEX/PCp/us-based  programs,  two  im¬ 
ages  were  used  from  the  existing  image  library.  The  j!7res4.img  wras  again  used,  this 
time  to  determine  the  accuracy  of  the  filter  programs.  Both  feat-id. c  and  sizeit.c  wrere 
operationally  tested  with  the  image  lOiu'grid.img  derived  from  a  standard  calibration 
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Figure  4.1 :  Unfiltered  jl7res4.img. 


image  [Ref.  2].  Resultant  data  acquired  from  these  images  was  compared  to  data  in 
Reference  4. 

Speckle  index  value  comparisons  weie  used  to  determine  the  effectiveness  of  the 
filter  programs.  In  each  case,  the  respective  filter  was  run  through  six  iterations  of 
the  full  screen  image,  jl7rcs4.img.  Table  4.4  contains  the  tabulation  of  the  results 
acquired,  and  the  results  from  the  original  C  programs  are  listed  in  Table  4.5.  A  file 
comparison  was  performed  to  ensure  that  the  image  used  for  testing  was  identical 
to  the  full  screen  image  used  previously.  Comparison  of  the  data  obtained  from  the 
current  iterations  to  data  from  Kaeser's  programs  indicates  slight  variations  in  the 
speckle  index  values.  The  initial  speckle  index  values  for  the  two  cases  differ.  One 
possible  factor  attributing  to  this  difference  is  the  inconsistency  of  integer  conversions 
between  the  C  language  programs  [Ref.  4].  The  more  significant  measure  of  filter 
performance  is  the  relative  decrease  in  speckle  index  upon  subsequent  filter  iterations. 
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TABLE  4.4:  Speckle  Index  Reduction  Results  Obtained  from  Modified 
Programs 


NUM 

2  SIGMA 

L  STAT. 

GEOMETRIC 

0 

0.293046 

0.293046 

0.293046 

1 

0.196282 

0.180061 

0.218982 

2 

0.144286 

0.125157 

0.174153 

3 

0.115754 

0.099063 

0.143760 

4 

0.098420 

0.084035 

0.121202 

5 

0.0870 68 

0.074063 

0.103756 

6 

0.079042 

0.066919 

0.089894 

TABLE  4.5:  Speckle  Index  Reduction  Results  (Ref.  4) 


NUM 

2  SIGMA 

L  STAT. 

GEOMETRIC 

0 

0.304005 

0.304005 

0.304005 

1 

0.153767 

0.194121 

0.224875 

2 

0.098967 

0.151951 

0.178222 

3 

0.074415 

0.129769 

0.146152 

4 

0.060615 

0.115938 

0.122347 

5 

0.051691 

0.106372 

0.103983 

6 

0.045354 

0.099367 

0.089425 

Therefore,  the  dissimilarities  of  initial  speckle  index  value  are  not  detrimental  to  filter 
function,  and  all  filter  programs  produce  desirable  reduction  of  speckle.  Consequently, 
all  evidence  indicates  that  the  ITEX/PCp/us-based  filter  programs  function  correctly. 
Figure  4.2  illustrates  jllrts^Amg  after  six  iterations  of  the  geometric  filter. 

The  results  obtained  from  the  ftatjd.c  and  sizeit.c  programs  were  identical  to 
those  obtained  from  Kaeser’s  programs  [Ref.  4].  The  SCALE-FACTOR  entry  in  the 
thesis. h  header  file  had  to  be  modified  in  order  for  the  output  files  from  the  respective 
sizing  programs  to  be  congruent.  Otherwise,  the  program  files  produced  the  same 
outputs  as  the  original  C  language  programs. 
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Figure  4.2:  jl7res4-img  Image  After  Six  Iterations  of  the  Geometric  Filter. 
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V.  CONCLUSIONS 


A.  SUMMARY  OF  FINDINGS 

This  thesis  accomplished  all  major  objectives  as  set  forth  in  the  introductory 
chapter.  No  major  problems  were  anticipated  with  the  integration  of  the  ITEX/P Cplus 
libraries  and  the  COMPAQ  386/20  computer  system.  Due  to  a  memory  conflict  cre¬ 
ated  within  the  COMPAQ,  the  PCVISIONp/us  had  to  be  installed  into  the  IBM 
PC/AT  to  complete  C  program  modifications.  Therefore,  the  bulk  of  this  research 
was  conducted  on  the  PC/AT  used  in  previous  work.  Although  several  of  the  pro¬ 
grams  produced  in  Kaeser’s  research  [Ref.  4]  were  modified  with  no  major  difficulties, 
type  definitions  created  major  faults  in  the  filter  routines  which  had  to  be  resolved. 
It  was  the  goal  of  this  thesis  to  improve  program  portability  and  optimization.  For 
the  most  part,  this  goal  was  achieved.  Executable  program  size  did  increase,  but  this 
was  the  result  of  program  structure  changes  and  library  modifications  introduced  by 
the  ITEX/Pcp/us  system  upgrade. 

B.  FUTURE  WORK 

Further  verification  testing  is  necessary  to  sufficiently  demonstrate  the  operation 
of  the  C  language  programs.  A  minimal  amount  of  experimentation  was  accomplished 
to  compare  the  acquired  results  with  data  obtained  from  the  previous  version  of  the 
programs.  Further  testing  is  mandatory  for  complete  program  validation  and  should 
be  pursued  in  subsequent  thesis  research  prior  to  porting  the  programs  over  to  other 
computer  systems.  Some  effort  was  expended  on  the  application  of  the  NDP  C-386 
compiler  to  eliminate  the  need  for  virtual  arrays  for  data  processing.  These  efforts 
were  largely  unsuccessful  due  to  the  incompatibility  of  the  ITEX/PCp/ws  libraries 
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with  the  subject  compiler  system.  This  area  of  research  could  be  further  investigated 
for  future  thesis  work.  The  following  topics  may  also  be  considered  for  future  work: 

•  Utilization  of  Macintosh  board  for  image  processing. 

•  Measurement  of  image  particles  utilizing  a  scanning  electron  microscope  for 
sizing  routine  validation. 

•  Further  validation  of  image  measurements  using  new  calibration  objects. 

•  Integration  of  SUN  workstations  for  image  processing. 


APPENDIX  A:  PROGRAM  HEADER  FILE: 

thesis. h 

/*  PURPOSE: 

To  completely  list  all  necessary  include  files,  manifest  constants, 
and  MACRO  definitions  required  by  all  the  programs  within  this 
package.  Any  changes  to  constants  in  thesis. h  will  guarantee 
redefinition  of  those  parameters  throughout  the  main  files. 


*/ 

/*  Include  files  for  use  with  ITEXPC  programs 
*/ 


# include 
•include 
* include 
•include 
•include 
•include 
•include 
•include 


<stdio.h> 
<stdlib.h> 
<math . h> 
<malloc.h> 
<gr aph . h> 
<conio.h> 
"itexpfg.h" 
"stdtyp.h" 


/*  Program  constants  used  in  the  C  program  nodules  of  the  ITEXPC 

Optical  System  used  for  Rocket  motor  Hologram  analysis.  •/ 

/*  Intial  ITEXPC  Board  Jumper  Settings  */ 

•define  MEHBASE  OxDOOOOL  /*  Iter  board  base  memory  start  address  */ 

•define  REGBASE  0x300  /*  Itex  board  base  register  start  address  */ 

•define  MEMORY  DUAL  /*  Itex  board  memory  type  */ 

/*  ITEXPC  initial  AOI  (area  of  interest)  settings  */ 


•define  IXS  0 

•define  IYS  0 

•define  IROW  480 

•define  ICOL  512 

•define  COLOR  150 

•define  XSIZE  512 

•define  YSIZE  512 

•define  DEPTH  8 


/*  ItexPC  LUT  (Look  Up  Table)  Variables  */ 
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/*  Threshold  Limits  */ 


•define 

LOWEST 

0 

/•  Equates  to  Black  for  lowcut  value  */ 

•def ine 

HIGHEST 

265 

/*  Equates  to  White  for  highcut  value  */ 

•define 

BUCK 

0 

•define 

WHITE 

255 

/*  Filter  array 

sizes  and 

limits  •/ 

•def ine 

SHUN 

25 

•def ine 

HUM 

9 

•define 

HIGH 

254 

•def ine 

LOW 

0 

/*  Sizing  magnification  factor  used  in  sizing  program  */ 

•define  SCALE.FACTOR  10.000000  /*  Based  on  (  317.6  /  85  )  */ 

/*  Virtual  array  Header  File  Definitions  •/ 

•define  header  7 


/*  Virtual  Array  Control  Block  typedef  */ 
typedef  struct  { 


FILE  efile; 
long  size; 
int  elsize; 
char  * buffer; 
int  buf_elsize; 
int  buf.size; 
char  *blank_rec; 


/*  pointer  to  file  control  block  */ 

/*  number  or  array  elements  in  file  */ 

/•  number  or  bytes  in  each  element  */ 

/*  pointer  to  array  buffer  */ 

/*  size  of  element  in  buffer  including  index  */ 
/•  number  of  elements  in  buffer  */ 

/*  pointer  to  initialization  record  */ 


} 

VACB 


/*  used  for  extending  file  */ 

;  /*  Virtual  Array  control  block  type  name  */ 


/*  Virtual  Array  Access  Prototypes  */ 

int  init_v_axray(char  *f ilename.int  rec_size,char  filchar); 
VACB  *open_v_array(char  *f ilename , int  buffer.size) ; 
void  close_v_array(VACB  *v_array) ; 
void  *access_v_rec(VACB  *v_array .long  index); 

/*  Virtual  Array  Access  Macros  */ 

•define  VREC(y)  ((items  *)access_v_rec(item_array ,y)) 
•define  gdesc(y)  VREC(y)->v_gdesc 
•define  G(x,y)  VREC(y)->v_gdesc[x] 

/*  Virtual  Array  element  structure  typedef  */ 
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typedef  struct  { 

unsigned  char  v.gdesc [ICOL] ; 

> 

iteas ; 
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APPENDIX  B:  PROGRAM  FILE:  geofil.c 


/•PURPOSE  :  Provides  lor  filtering  an  inage  nsing  the  geometric  filter 
algorithm.  This  program  processes  the  image  and  provides 
for  the  operator  to  select  the  number  of  iterations. 

The  area  of  filtering  for  the  image  is  controlled  by  the 
operator  by  keyboard  input  of  number  of  ROWS  and  COLS  to 
process.  Program  sill  alios  for  a  Maximum  input  of  480  ross 
and  512  cols  sith  a  minimum  of  1  ros  and  1  column. 

The  resulting  filtered  image  pixel  value  is  stored  on 
disk  in  a  virtual  array,  that  is  sritten  to  the  image 
processing  screen  shen  all  calculations  are  complete. 

Total  calculation  time  for  a  512x480  image  is  3403.13  sec. 

DISK  SPACE  REQUIRED  for  virtual  array  storage  is  250  KB  minimum 
available  on  the  disk  at  start  of  the  routine. 


*/ 

kinclude  "thesis. h" 

char  filename [20] ,  comline [200] ; 

int  sros,  scol,  times,  ros,  col; 

main() 

{ 

int  flag; 

extern  int  ros,  col,  sros,  scol,  times; 

printf ("\n\tReady  to  Load  IMAGE  FILE?. . .Yes  (y)  /  Ho  (n)  "); 
flag=getch() ; 

if (flag  ==  ' Y ’ II  flag  ==  ’y') 

{ 

startitO ; 
readitO ; 

> 

if(  flag  ==  ’H’  ||  flag  ==  'n’  )  return(O); 

geti("\n\nEnter  number  of  ross  (Max=480)  to  use  in  SPECKLE  CALC  :  ".fcsros); 
geti ( "\nEnter  number  of  columns  (Max=512)  to  us  in  SPECKLE  CALC  :  ",Ascol); 
geti("\n\n\nEnter  number  of  filter  iterations  to  do. . . " .Atimes) ; 
geti("\nEnter  the  number  of  IMAGE  ross  (MAX  =  480)  to  FILTER  :  ",Aros); 
geti("\nEnter  the  number  of  IMAGE  cols  (MAX  =  512)  to  FILTER  :  M,tcol); 
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switch(times) 

{ 

case  0: 

apeckleO ; 

printf ("\nSave  image  to  Diak  File?.. Yea  (y)  /  lo  (n)  "); 
flag=getch() ; 

if (flag  ==  *Y'II  flag  ==  *y*)  aaveitO; 

ezit(l) ; 

break; 

caae  1 : 

break; 

default : 

break; 

> 

do 

{ 

apeckleO ; 
geofilO ; 

> 

while(times  ==  0); 

printf (”\nSave  image  to  Diak  File?.. Yea  (y)  /  lo  (n)  "); 
flag=getch() ; 

if  (flag  ==  *  Y  *  I  I  flag  ==  ’y’)  aaveitO; 


/*  geometric  filter  algorithm  */ 

geofilO 

i 

VACB  eitem_ array ; 

extern  int  row,  col,  tinea; 

int  fl,  f2,  gl.  g2 ,  g3 ,  x,  y; 

int  a,  b,  c,  d,  pixel,  calla; 

/*  create  a  virtual  array  for  the  array  of  filtered 
pixel  values  the  size  of  the  image  (  480x512  )  */ 

init_v_array("ITEMS. VAR" .sizeof (items) ,IULL) ; 

/*  open  the  virtual  array,  reserve  buffer  apace  for  100  elements  */ 
item.array  =  open_v_array("ITEMS. VAR", 100) ; 

/*  start  local  statistic  filter  routine  */ 


if (  row  >  BROW  )  row  =  BROW; 
if (  col  >  BC0L  )  col  =  BC0L; 


for(  y  =  1  ;  y  <  row  ;  y++  )  /•  Zero  the  initial  array  */ 

{ 

for(  x  =  1  ;  x  <  col  ;  x++  ) 

{ 

G(((long)x) , ((long)y))  =  IULL; 

> 

> 

for (calls  ~  tines;  calls  !=  0;  calls — ) 
c  *  4; 

printf ( "\n\n\tF iltsr  Running  — >  y,d  Runs  aftsr  this  ons! . . .\nM,calls-l) ; 
ior(  ;  c  ;  ) 
i 

svitch(c) 

{ 

case  4  : 

a=l ;  b=0;  c=3;  d=l; 
printf ("\nc  =  3  "); 

tl(a,b,c,d,item_array) ; 
break ; 

case  3: 

a=0 ; 
b=l ; 
c=2; 

printf ("\nc  =  3  "); 

tl(a,b,c,d,item_array) ; 
break; 

case  2: 

a=l ; 
b=l ; 
c=l; 

printf ("\nc  =  2  ”) ; 

tl(a,b,c,d,iteiB_array) ; 
break; 

case  1 : 

a=l; 
b=-l ; 
c=0; 

printf ("\nc  =  1  "); 

tl(a,b,c,d,itere_array) ; 
printf ("\nc  =  0  EHD") ; 
speckle()  ; 
brer*- 


default  :  break; 

> 

} 

> 

/*  close  the  virtual  array  */ 
elose_v_array(iten_array) ; 


/•  subroutine  for  geofil  for  steps  1  and  2  •/ 


tl(a,b,c,d,iteB_array) 
int  a,  b,  c,  d; 
VACB  *iten_ array; 

{ 


int 

int 

long 

extern  int 


fl.  f 2 ,  gl.  g2 ,  g3; 
pixel,  x,  y; 

*1.  yi; 

row, col, tines; 


for(  y  =  1  ;  y  <  row-1  ;  y+>  ) 

{ 

for(  x  =  1  ;  x  <  col-1  ;  x++  ) 

{ 

xl  =  ((long)  x); 
yl  =  ((long)  y); 
fl  =  rpixel(  x,y  ); 
f 2  =  rpixel(  x-a,y-b  ); 

G(xl,yl)=  nax(  fl,  nin(  f2-l,  fl+1  )  ); 

> 

> 

printf("  step  1"); 

for(  y  =  1  ;  y  <  row-1  ;  y++  ) 

{ 

for(  x  =  1  ;  x  <  col-1  ;  x++  ) 

{ 

xl  =  ((long)  x); 
yl  =  ((long)  y) ; 
gl  =  G(  xl.yl); 

g2  =  G(((long) (x-a) ) , ((long) (y-b))) ; 

g3  =  G(( (long) (x+a)), ((long) (y+b))); 

pixel  =  max(  gl,  nin(  g2,  nin(  gl+1,  g3+l  )  )  ) 

wpixel(x,y .pixel) ; 

} 

> 

printfC  step  2"); 
if (  d==l  ) 

{ 

a=  -a; 
b=  -b; 
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d=  0; 

tl(a,b,c,d,itea_array) ; 


> 

else  if(  d==0  ) 
d=l; 

t2(a,b,c,d,iteB_ array) ; 
return ; 

> 

> 


/*  subroutine  1 or  geolil  lor  steps  3  and  4  */ 
t2(a,b,c,d,item_array) 
int  a,  b,  c,  d; 

VACB  *item_array; 


< 


int 

int 

long 

extern  int 


11.  12.  gl.  g2,  g3; 
x,  y.  pixel; 
xl.  yi; 

row,  col,  times; 


lor(  y  =  1  ;  y  <  row-1  ;  y++  ) 

{ 

lor(  x  =  1  ;  x  <  col-1  ;  x++  ) 

{ 

xl  =  ((long)  x); 

yl  =  ((long)  y); 

11  =  rpixel(  x.y  ); 

12  =  rpixel(  x-a,y-b  ); 

G(xl ,yl)=  min(  11,  max(  12-1,  11+1  )  ); 

> 

> 

print! ("  step  3"); 

1 or(  y  =  1  ;  y  <  row-1  ;  y++  ) 

{ 

lor (  x  =  1  ;  x  <  col-1  ;  x++  ) 

xl  =  ((long)  x); 
yl  =  ((long)  y); 
gl  =  G(  xl.yl); 

g2  =  G(  ( (long) (x-a) ) , ((long) (y-b) ) ) ; 

g3  =  G(((long) (x+a)) , ((long) (y+b))) ; 

pixel  =  max(  gl.  min(  g2,  min(  gl+1,  g3+l  )  )  ); 

wpixel(x,y .pixel) ; 

> 

> 

print! ("  step  4") ; 
il(  d==l  ) 


a=  -a; 
b=  -b; 
d=  0; 

t2(a,b,c,dl iten.array) ; 

> 

•Isa  i t(  d==0  ) 

{ 

d  =  1; 
return; 

> 

> 
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APPENDIX  C:  PROGRAM  FILE:  lstat.c 


/♦PURPOSE  :  Provides  lor  filtering  an  iaage  using  the  local  statistical 

algorithm.  This  prograa  processes  the  iaage  by  using  a  local 
SzS  array  of  pixels  to  calculate  a  statistical  value  for  the 
central  pixel  of  the  local  array. 

The  area  of  filtering  for  the  iaage  is  controlled  by  the 
operator  by  keyboard  input  of  nuaber  of  ROWS  and  COLS  to 
process.  Prograa  sill  alios  for  a  Haxiaua  input  of  480  roes 
and  S12  cols  sith  a  ainiaua  of  1  ros  and  1  column.  Prograa 
also  requires  the  iaage  Std.  Dev.  as  calculated  by  the  genfunc.c 
Program.  The  resulting  filtered  image  pixel  value  is  stored  on 
disk  in  a  virtual  array,  that  is  written  to  the  iaage  processing 
screen  when  all  calculations  are  complete. 

Total  calculation  time  for  a  512x480  iaage  is  1324.47  sec. 

DISK  SPACE  REQUIRED  for  virtual  array  storage  is  250  FB 
minimum  available  on  the  disk  at  the  start  of  the  routine. 

*/ 

kinclude  "thesis. h" 

char  filename [20] ,  comline [200] ; 

int  row,  col,  sros,  scol,  times; 

float  devl; 

main() 

{ 

int  flag; 

extern  int  ros,  col,  times; 

pr intf ( "\n\tReady  to  Load  IMAGE  FILE?... Yes  (y)  /  lo  (n)  "); 
f lag=getch() ; 

if (flag  ==  ’ Y ’  I  I  flag  ==  ’y') 

< 

startitO  ; 
readitO ; 

> 

if(flag  ==  ’I’  II  flag  ==  ’n’)  return(O); 

geti("\n\nEnter  number  of  ross  (Max=480)  to  use  in  SPECKLE  CALC  :  ",Asros); 
geti("\nEnter  number  of  columns  (Max=512)  to  us  in  SPECKLE  CALC  :  ",Ascol); 
geti("\n\n\nEnter  number  of  filter  iterations  to  do. . . " .Atimes) ; 


geti("\nEnter  the  number  of  IMAGE  rose  (MAX  =  480)  to  FILTER  :  ",Arow); 
geti("\nEnter  the  number  of  IMAGE  cols  (MAX  =  612)  to  FILTER  :  ",Acol); 
dev(Adevl) ; 

printf ("\nThe  standard  deviation  for  this  image  is  %.2f\n" ,devl) ; 


do 

{ 

■peckleO ; 

IstatO ; 

> 

while(tines  ==  0); 

printf("\nSave  inage  to  Disk  File?.. Yes  (y)  /  lo  (n)  "); 
flag-getch() ; 

if  (flag  ==  *  Y '  1 1  flag  ==  'y')  saveitQ; 


/* . */ 

IstatO 


VACB 

long 

extern  float 

extern  int 

register  int 

int 

float 

float 

int 


*iten_array ; 

xi.yi; 

devl ; 

row,  col,  tines,  srow,  scol; 
n,  n,  nn; 

pixel, a, b.x.y, calls; 

suml,  sum2,  svar,  lvar,  tvar; 

lnean,  lnean2,  stddev,  stddev2,  k; 

IdataCSIUM];  /*  SIUM  defined  in  Header  file  */ 


/*  create  a  virtual  array  for  the  array  of  filtered 
pixel  values  the  size  of  the  inage  (  480x512  )  */ 

init_v_array(" ITEMS. VAR",sizeof (itens) ,IULL) ; 

/*  open  the  virtual  array,  reserve  buffer  space  for  100  elenents  */ 
item.array  =  open_v_array("ITEMS.VAR", 100) ; 

/*  start  local  statistic  filter  routine  */ 


if (  row  >  IR0V  )  row  =  IR0W; 
if (  col  >  I COL  )  col  =  XC0L; 

stddev  =  devl  /  2SS.0  ;  /*  lomalize  deviation  value  */ 

stddev2  =  stddev  *  stddev  ; 

for(calls  =  tines;  calls  !=  0;  calls — )  { 

printf ("\n\n\tFilter  Running  — >  %d  Runs  after  this  one! . . ,\n",calls-l) ; 
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for  (  y  =  2  ;  y  <  row  -  2  ;  y++) 

{ 

for  (  x  =  2  ;  x  <  col  -  2  ;  x++) 

{ 

pixel  =  rpixel(  x,y  ); 
sunl  =  0; 
sum2  *  0; 
nn  =  0; 

for(  n  =  y-  2;n<y  +  3;  n++  ) 

{ 

for(  n=x-2;n<x  +  3;  »+♦  ) 

{ 

ldata[nn]  =  rpixel(  a,n  ); 
surnl  ♦=  IdataCnn]; 
nn++; 

> 

> 

lnean  =  sunl  /  SIUM  ; 

lnean2  =  lnean  *  lnean  ; 

for(  nn  =  0  ;  nn  <  SIUM  ;  nn++  ) 

{ 

avar  =  (IdataCnn]  -  lnean  )  *  (  IdataCnn]  -  lnean  ); 
eun2  +=  svar ; 

> 

lvar  =  sum 2  /  SIUM  ; 

tvar  =  fabs(  (  (lvar  +  lnean2)  /  (stddev2  ♦  1)  )  -  lmean2  ) 
k  =  tvar  /  (  (  stddev2  *  lnean2  )  ♦  tvar  )  ; 
xl  =  ((long)  x); 
yl  *  ((long)  y); 

G(  xl,yl  )  =  (int)  (  lnean  ♦  k  *  (  pixel  -  lnean  )  )  ; 

> 

> 

/*  write  contents  of  the  G  array  to  image  screen  */ 

printf ("\n\a\tWriting  filtered  inage  to  screen..."); 

for(  b  =  2  ;  b  <  row  -  2  ;  b++  ) 

( 

for(  a  =  2  ;  a  <  col  -  2  ;  a++  ) 

wpixel(  a,b,G(((long)a) , ( (long)b) ) ) ; 

> 

} 

speckle() ; 

/*  close  the  virtual  array  */ 
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cloaa_v_array(itaa_array) ; 


APPENDIX  D:  PROGRAM  FILE:  2sigma.c 


/♦PURPOSE  :  Provides  f or  filtering  an  iaage  using  a  2-sigma  statistical 

algorithm.  This  program  processes  the  in age  by  using  a  local 
SzS  array  of  pixels  to  calculate  a  statistical  value  for  the 
central  pixel  of  the  local  array.  The  algoritha  vill  eliminate 
pixels  froa  the  suanation  that  are  greater  than  2-Sigaa  in 
value  froa  that  of  the  central  pixel. 

The  area  of  filtering  for  the  iaage  is  controlled  by  the 
operator  by  keyboard  input  of  nuaber  of  ROWS  and  COLS  to 
process.  Program  will  allow  for  a  Maximum  input  of  480  rows 
and  512  cols  with  a  Minimum  of  1  row  and  1  column.  Program 
also  requires  the  iaage  Standard  deviation  as  obtained  using 
genfunc.c  program.  The  resulting  filtered  pixel  value  is 
stored  on  disk  in  a  virtual  array  and  written  to  screen  later. 

Total  calculation  tine  for  a  512x480  iaage  is  760.55  sec 

DISK  SPACE  REQUIRED  for  virtual  array  storage  is  250  KB 
minimum  available  on  the  disk  at  start  of  the  routine. 


♦/ 


•include 

"thesis .h" 

char 

f ilename(20] ,  comline [200]  ; 

int 

row,  col.t  imes,  srow,  scol 

float 

devl ; 

main() 

int 

flag; 

extern  int  row,  col,  times; 

printf ("\n\tReady  to  Load  IMAGE  FILE  from  disk?... Yes  (y)  /  Mo  (n)  "); 
f lag=getch( ) ; 

if (flag  ==  * Y *  I  I  flag  ==  'y’) 

{ 

startitO  ; 
readit () ; 

> 

if(flag  ==  I!  flag  ==  'n')  return(O); 

geti("\n\nEnter  number  of  rows  (Max=480)  to  use  in  SPECKLE  CALC  :  ".tsrow); 
geti("\nEnter  nuaber  of  columns  (Kax=5l2)  to  us  in  SPECKLE  CALC  :  " .tscol) 
geti("\n\n\nEnter  number  of  filter  iterations  to  do. . . " .Atimes) ; 


geti("\nEnter  the  number  of  IMAGE  rows  (MAX  =  480)  to  FILTER  :  ",Arow); 
geti("\nEat#r  the  number  of  IMAGE  cols  (MAX  -  512)  to  FILTER  :  ",*col); 
dev(ftdevl); 

printf ("\n\tThe  standard  deviation  for  this  image  is:  %.2f\n" ,devl) ; 
do 

{. 

speckleO; 
signaQ ; 

> 

while(  tines  ==  0  ) ; 

printf ("\nSave  inage  to  Disk  File?.. Yes  (y)  /  lo  (n)  "); 
flag-getchO; 

if  (flag  ==  *  Y  ’  I  I  flag  ==  Jy*)  saveitO; 


> 

Zee*****************************************************************/ 


sigmaO 

{ 

VACB 

extern  int 
extern  float 
long 

register  int 

int 

int 

float 

int 


*item_array ; 

row,  col,  tines,  srow,  scol; 
devl; 

Xi.  yi; 

n,  n,  nn; 

pixel,  sunl,  sun2,  delta,  high; 
low,  a,  b,  x,  y,  calls; 
hvar,  lvar,  stddev; 

IdatatSIUM] ;  /*  SWUM  defined  in  Header  file  */ 


/*  create  a  virtual  array  for  the  array  of  filtered 

pixel  values  the  size  of  the  image  (  512x480  )  */ 


init_v_array( “ITEMS. VAR", sizeof (items) ,IULL) ; 


/*  open  the  virtual  array,  reserve  buffer  space  for  100  elements  */ 
item_array  =  open_v_array("ITEMS.VAR",100) ; 


/*  start  local  statistic  filter  routine  */ 


if (  row  >  *R0W  )  row  =  IROV; 
if (  col  >  ICOL  )  col  =  *C0L; 

stddev  =  devl  /  255.0  ;  /*  formalize  deviation  value  */ 

hvar  =  1.  +2.  *  stddev  ;  /*  find  high  and  low  sigma  */ 

lvar  =  1.  -  2.  *  stddev  ; 

for  (calls  =  times;  calls  !=  0;  calls — )  { 

printf  ("\n\n\tFilter  running  -->  */,d  Runs  after  this  one  !  .  .  .  \n"  ,calls-l) ; 
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for  (  y  =  2  ;  y  <  row  -  2  ;  y++) 

•C 

lor  (  x  =  2  ;  x  <  col  -  2  ;  x++) 

{ 

pixel  =  rpixel(  x,y  ); 

s uni  =  0; 

sun2  =  0; 

an  =  0; 

delta  =  0; 

high  =  (int)(hvar  *  pixel); 
low  =  (int)(lvar  *  pixel); 

for(  n=y-2;n<y  +  3;  n++  ) 

lor(  m=x-2;n<x  +  3;  m++  ) 

{ 

xl  =  ((long)  x); 
yl  =  ((long)  y); 

IdataCnn]  =  rpixel(  m,n  ); 

il(  (  low  <=  IdataCnn]  )  kk  (  IdataCnn]  <=  high  )  ) 

{ 

s uni  +=  IdataCnn] ; 
delta++; 

> 

nn++ ; 

> 

> 

iKdelta  <=  2)  /*  correct  shot  noise  —  4  neighbor  average  */ 

{ 

sun2  =  (  rpixel(x,y-l)+rpixel(x,y+l)+rpixel(x-l ,y)+rpixel(x+l ,y)  ); 

G(  xl.yl  )  =  (pixel  +  s\im2)  /  5; 

continue; 

> 

G(  xl.yl  )  =  suml  /  delta  ; 

> 

> 

print! ("\n\tWriting  filtered  inage  to  screen...’*); 

for(  b  =  2  ;  b  <  row  -  2  ;  b++  ) 

{ 

for(  a  =  2  ;  a  <  col  -  2  ;  a++  ) 

wpixel(  a,b,G(((long)a) ,((long)b))  ); 

> 

speckleO  ; 


/*  close  the  virtual  array  */ 


clos«_ v .array ( itaa.array ) 


APPENDIX  E:  PROGRAM  FILE:  threshit.c 


/♦PURPOSE  :  This  program  thresholds  the  image  on  screen  by  taking 
an  operator  input  value  and  forcing  all  image  pixel 
values  above  the  threshold  value  to  the  value  of  BLACK 
and  all  those  belos  the  value  to  the  WHITE  value. 

(  Background  ==  WHITE  ;  Feature  ==  BLACK  ) 


*/ 

•include  "thesis. h" 

char  filename [20] ,  coaline[200] ; 

int  sros,  scol; 

mainO 

{ 

int  flag; 

printf ("\nThis  program  sill  threshold  the  input  image  desired."); 
printf ("\n\n\tReady  to  Load  Image?... Yes  (y)  /  lo  (n)  "); 
flag=getch(); 

if (flag  ==  *Y*  I  I  flag  «  *y* ) 

startitO  ; 
readitQ  ; 

> 

if (flag  ==  >I*| I  flag  ==  >n>) 

initializeO ; 
return(O) ; 

> 

printf ("\n\n\tTHRESH0LD  the  image?... Yes  (y)  /  Ho  (n)  "); 
flag  =  getch(); 

if  (flag  ==  >Y>  I  I  flag  ==  >y>)  threshitO; 


> 
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APPENDIX  F:  PROGRAM  HEADER  FILE 

feat  Jd.c 


/♦PURPOSE  :  Labels  and  identifies  each  feature  in  an  image  on  the  ITEX 
system.  Reads  pixel-by-pixel  and  groups  the  objects  by 
assigning  a  unique  ID  number  to  each  feature  or  object  so 
that  they  can  be  processed  by  other  routines. 

The  ID  routine  requires  a  thresholded  on  image  screen  image 
or  input  of  a  previously  saved  thresholded  image  from  disk  or 
imput  from  disk  and  then  thresholding.  This  program  module 
includes  an  optional  call  to  threshold  if  desired. 

(  Background  ==  WHITE  ;  Feature  ==  BLACK  ) 

*/ 

finclude  "thesis. h" 

char  filename [20] ,  comline[200] ; 

int  srow,  scol; 

mainO 

< 

int  flag; 

printf ("\n\aYou  MUST  USE  A  THRESHOLDED  IMAGE  for  this  Program!!"); 
printf ("\n\n\tReady  to  Load  image  from  disk?... Yes  (y)  /  lo  (n)  "); 
flag=getch(); 

if (flag  ==  *Y» I  I  flag  ==  >y’) 

{ 

startitO  ; 
readitO ; 

} 

if  (flag  ==  'I*  II  flag  ==  'n'  )  retum(O); 

printf ("\n\n\t*eed  to  THRESHOLD  the  image?... Yes  (y)  /  lo  (n)  "); 
flag  =  getchO; 

if  (flag  ==  ’Y’ll  flag  ==  ’y’)  threshO; 
f eat_id() ; 

printf ("\n\nSave  image  to  Disk  File?.. Yes  (y)  /  to  (n)  "); 
flag=getch() ; 

if  (flag  ==  'Y’ll  flag  ==  'y')  saveitQ; 

> 


/* 

inage  feature 

identification  algorithm  */ 

feat 

_id() 

register  int 

x.  y; 

int 

a,  xl.  yl.  x2,  y2.  x3.  y3,  x4.  y4; 

int 

fid  =  0; 

int 

gid  =  0; 

int 

teap.fid  =  0; 

int 

te*p_gid  =  0; 

int 

staxfl,  nl,  n2,  nla,  nib.  nlc,  kk,  kkl,  nnl; 

print!  ("\n\n\tSTEP  OIE  in  Progress  ...**)  ; 

ior(  y  =  1  ;  y  <  *ROW  ;  y++  ) 

{ 

lor(  x  =  0  ;  x  <  ICOL  ;  x++  ) 

{ 

if(  rpixel(x.y)  ==  WHITE  }  continue; 
if  (  x  ==  0  ) 

{ 

wpixel(  x.y.fid  ); 
continue ; 

> 

il(  (  rpixel(  x-l.y  )  !=  WHITE  )  ) 

vpixel(  x,y,rpixel(x-l,y)  ); 
continue; 

> 

if (  rpixel(x.y-l)  !=  WHITE  ) 

{ 

upixelC  x,y,rpixel(  x,y-l  )  ); 
continue; 

> 

wpixel(  x.y.fid  ); 
a  =  x  +  1 ; 

while(l) 

if(  rpixel(  a.y  )  ==  WHITE  ) 

{ 

wpixel(  x.y.fid  ); 
step(  kfid.kgid  ); 
break; 

} 

if(  rpixel(a.y-l)  •=  WHITE  ) 

{ 

wpixel(  x,y,rpixel(  a,y-l  )  ); 
break ; 

> 
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a  +=  1; 


> 

>  /*  and  x  -  for  */ 

>  /*  and  y  -  for  */ 

fid  ♦=  gid*HIGH; 

printf ("\n\n\tSTART  STEP  TWO...  "); 

printf ("\n\nEntar  tha  BEST  estimate  of  tha  aax  feature  length"); 
geti("\n\nlt'8  Battar  to  ba  too  large  than  too  small !!" .kmaxfl) ; 
printf ("\n\n\tRUIIIIG") ; 

for(  y  =  1  ;  y  <  IROW  ;  y++  ) 

for(  x  =  1  ;  x  <  ICOL  ;  x++  ) 

nl  =  rpixel(  x,y  ) ; 
n2  =  rpixel(  x,y-l  ); 

if (  (  nl  II  n2  ==  WHITE  )  II  (nl  ==  n2  ))  continue; 

x3  =  x  -  maxfl; 
x4  =  x  +  maxfl; 
y3  =  y  -  maxfl; 
y4  =  y  +  maxfl; 

if (  x3  <  LOW)  x3  =  LOW; 
if (  x4  >  ICOL  )  x4  =  ICOL; 

if (  y3  <  LOW)  y3  =  LOW; 

if (  y4  >  IROW  )  y4  =  IROW; 

for  (  yl  =  y3  ;  yl  <  y4  ;  yl++  ) 

< 

for  (  xl  =  x3  ;  xl  <  x4  ;  xl++) 

{ 

nla  =  rpixelC  xl ,yl  ) ; 

if  (  nla  ==  WHITE)  continue; 

if  (  nla  ==  nl  )  wpixel(  xl,yl,n2  ); 

> 

fid  -=  1; 

> 

> 

> 

printf ("\n\n\tSTEP  THREE  in  Progress"); 

for  (  y  =  1  ;  y  <  IROW  ;  y++  ) 

for(  x  =  1  ;  x  <  ICOL  ;  x++  ) 

{ 

kk  =  rpixel(  x,y  ) ; 

if(  kk  ==  WHITE  )  continue; 


56 


kkl  =  kk; 

nnl  =  temp_fid  +  temp.gid  *  HIGH; 
if(  kkl  <  nnl  )  continue; 
if (  (kkl  -  nnl)  <  -350  )  kkl  =  kk  +  508; 
if (  (kkl  -  nnl)  >=  -350  kk  (kkl  -  nnl)  <  -160  ) 
kkl  =  kk  +  254; 
if (  kkl  ! =  nnl  ) 

{ 

x3  =  x  -  maxfl; 

x4  =  x  +  maxfl; 

y4  =  y  +  maxfl; 

if (  x3  <  0)  x3  =  LOW; 

if (  x4  >  ICQL-1  )  x4  =  ICOL  -1; 

if(  y4  >  IROW-1  )  y4  =  IROW  -1; 

for(  yl  =  y  ;  yl  <  y4  ;  yl++  ) 

{ 

for(  xl  =  x3  ;  xl  <  x4  ;  xl++) 

{ 

nlc  =  rpixel(  xl.yl  ); 

if  (  (  nlc  ==  WHITE  )  I  I  (  nlc  !=  kk  )  Continue 
wpixel(  xl ,yl ,temp_f id  ); 

> 

> 

step(  ktemp_f id , kt emp_gid  ); 

> 

> 

> 

printf ("\n\n\t\tFEATURE  COUNT  IS:  ^".fid); 

maplut (  INPUT,  0,  IXS,  IYS,  NCOL,  IROW  ); 


APPENDIX  G:  PROGRAM  FILE:  sizeit.c 


/♦PURPOSE  This  routine  uses  an  existing  or  saved  image  that  has  been 

thresholded  and  then  ID’d  with  the  ID  Algorithm.  The  output 
of  this  program  is  a  tabular  output  of  the  calculated  Area, 
X-chord,  and  Y-chord  that  is  suitable  for  input  to  a 
Statistical  Analysis  Program.  Output  file  is  written  to 
default  drive  as  SIZE.DAT. 

C  Background  ==  WHITE  ;  Feature  ==  BLACK  ) 

♦/ 

^include  "thesis. h" 

char  filename [20] ,  comline[200]  ; 

int  srow,  scol; 

main() 

{ 

int  flag; 

printf ("\n\aUSE  IMAGE  processed  by  the  FEATURE  ID  program!!"); 

printf ("\n\tReady  to  Load  ID’d  image  from  disk?... Yes  (y)  /  Ho  (n)  "); 
flag=getch() ; 

if (flag  ==  *Y’ I  I  flag  ==  ’y’) 

{ 

start it () ; 
readitQ ; 

> 

if(  flag  ==  ’H’  II  flag  ==  'n'  )  return(O); 

initializeQ ; 
size_it() ; 

printf ("\n\tSave  image  to  Disk  File?... Yes  (y)  /  Mo  (n)  "); 

flag=getch() ; 

if  (flag  ==  'Y’ll  flag  ==  ’y’)  saveitO; 

> 

/♦  image  feature  sizing  algorithm  ♦/ 

size_it() 

{ 

register  int  nn,  x,  y,  xa,  ya,  xb,  yb; 

int  j,  ip,  ipa,  ipb,  xl,  yl,  x2,  y2,  maxfl,  num; 
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int 

xmax,  ymax,  xmin,  ymin,  maxarea,  minarea; 

int 

xflag=0; 

int 

yflag=0; 

int 

aflag=0; 

float 

area,  xchord,  ychord,  factor; 

float 

fxmax,  f ymax ,  fxmin,  fymin,  f maxarea,  fminarea 

float 

ci_l,  cf_2,  mag; 

long 

♦aptr ; 

long 

♦xptr; 

long 

♦yptr; 

long 

♦stop; 

FILE 

♦featdata; 

geti("\nEnter  the  number  of  Features" .fcnurn) ; 
getf ("\nEnter  Magnification  Factor" ,kmag); 
geti("\nEnter  the  Maximum  Feature  Length" ,ftmaxfl) ; 

/♦  Set  initial  test  nin/max  values  */ 

xmax  =  0 ; 

ymax  =  0 ; 

naxarea  =  0; 
xmin  *  512; 

ymin  =  512; 

minarea  =  512; 

/♦  Dynamically  allocate  memory  space  for  the  number  of 
features  entered  by  the  operator. 

Note:  This  eliminates  the  need  to  declare  large  arrays 
for  these  values  and  allows  maximum  number  of  features 

xptr  =  (  long  ♦  )  calloc(  num  .sizeof (long)) ; 

yptr  =  (  long  *)  calloc(  num  .sizeof (long)) ; 

aptr  =  (  long  *)  calloc(  num  .sizeof (long)) ; 

/*  Chech  for  successfull  memory  allocation  */ 

if (  ! xptr  1 1  ! yptr  1 1  ! aptr  ) 

{ 

printf ("\n0ut  of  Memory!!!  \n"); 

printf ("\nYou  may  have  entered  too  many  features  An") ; 
printf ("\nCorrect  problem  and  try  again."); 
exit(0) ; 

> 

/•  Calculate  Magnification  Constants  */ 

scale(kfactor) ; 

cf _1  =  (  factor  /  mag  ) ; 

cf_2  =  (  cf_l  *  cf_l  ); 


/*  Above  Conversion  factor  converts  size  to  microns 

and  is  related  to  the  input  magnification  for  the 
hologram.  The  defined  "  SCALE.FACTOR  "  constant 
value  must  be  set  for  the  proper  value  based  on  the 
equipment  set  up  during  actual  image  aquisition. 

The  constant  term  SCALE.FACTOR  is  defined  in  the 
Header  File  THESIS. H 

*/ 

/*  Begin  Sizing  Routine  */ 

printf ("\n\n\tStarting  to  SIZE!...*'); 

nn  =  0;  /*  Feature  number  counter  */ 

for(  y  =  0  ;  y  <  MROV  ;  y++  ) 

{ 

if (  nn  ==  num  ) break ;  /*  Quit  ehen  all  features  sized  */ 

for(  x  =  0  ;  x  <  HCOL  ;  x++  ) 

{ 

ip  =  rpixel(x.y); 

if (  (ip  ==  WHITE)  II  (ip  ==  BLACK)  Continue; 

if (  nn  ==  num  )break;  /*  Quit  when  all  features  sized  */ 

xl  =  x  -  maxfl;  /*  Set  up  sizing  Box  */ 

yi  =  y; 

x2  *  x  +  maxfl; 

y2  =  y  +  maxfl; 

if (  xl  <  LOW)  xl  =  LOW; 
if (  yl  <  LOW)  yl  =  LOW; 
if (  x2  >  ICOL  )  x2  =  ICOL; 

if (  y2  >  HROW  )  y2  =  IROW; 

for(  ya  =  yl  ;  ya  <  y2  ;  ya++  ) 

{ 

for(  xa  =  xl  ;  xa  <  x2  ;  xa++  ) 

{ 

ipa  =  rpixel(xa.ya) ; 
if(  ipa  ! =  ip  ) continue ; 
aflag++; 
xflag++; 

> 

if(  xf lag  >  *(xptr  +  nn  )  )  *(xptr  +  nn)  =  xflag; 
xflag  =  0; 

> 

*(aptr  +  nn)  =  aflag; 
aflag  =  0; 
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for(  xb  =  xl  ;  xb  <  x2  ;  xb++) 

{ 

f or(  yb  =  yl  ;  yb  <  y2  ;  yb++  ) 

ipb  =  rpixel(xb.yb) ; 

if(  ipb  !=  ip  )  continue; 

yflag++; 

wpixel(xb,yb, BLACK) ; 

> 

if(  yf lag  >  *(yptr  +  nn  )  )  *(yptr  +  nn)  =  yflag; 
yflag  =  0; 

> 

/*  Calculate  Min/Max  values  */ 

xmax  s  max(  xnax,  * (xptr  +  nn)  ); 

ymax  =  max(  ymax,  *(yptr  +  nn)  ); 

max area  =  max (  maxarea,  *(aptr  ♦  nn)  )  ; 
xmin  =  min(  xmin,  *(xptr  +  nn)  )  ; 

ymin  =  min(  ymin,  *(yptr  +  nn)  )  ; 

m inarea  =  min(  ainarea,  *(aptr  +  nn)  )  ; 

nn++;  /*  increment  counter  */ 

printf  ("\n\n\tFeature  */,d  complete. .  .M,nn); 

> 

>  /*....  end  initial  x/y  ior  loop  ....*/ 

/+  output  data  section  */ 

printf ("\n\a  Sending  output  data  to  screen  and  FILE  SIZE.DAT  !"); 

/*  vaitO ;  ♦/ 

if(  (  featdata  =  fopen("size.dat","w")  )  ==  MULL  ) 

{ 

printf ("CAMMOT  Open  output  file  '/,s\n",  efeatdata); 
return(l) ; 

} 

printf ("  ID  MO  AREA  X-Width  Y-Midth  \n"); 

for(  j  =  0  ;  j  <  num  ;  j++  ) 

{ 

fprintf(  featdata, '"/,10d  */.10.3f  '/,10.3f  */.10.3f\n",  j+1,  *(aptr  +  j)  *  cf_2, 
*(xptr  +  j)  *  cf_l,  *(yptr  +  j)  *  cf_l  ); 

printf ("\ny.lOd  V.10.3f  ’/,10.3f  */,10.3f  ",  j  +  1,  *(aptr  +  j)  *  cf_2, 

*(xptr  +  j)  *  cf_l,  *(yptr  +  j)  *  cf_l  ); 

} 

/*  This  section  can  be  used  to  put  total  information  on  bottom  of  data  file 
if  desired.  File  written  now  to  be  used  with  STATGRAPHICS  */ 


61 


/*  fprintf deatdata, "\n\n\nMax  X-chord=  Xt  Max  Y-chord=  'It  Max  Area=  'it  ", 
x»ax*cf_l ,y*ax*cf_l ,aaxarea*cf _2) ; 

fprintf deatdata, "\nMin  X-chord=  Xt  Min  Y-chord=  Xi  Min  Area=  'it  ", 
xain*cf_l  ,jr»in*cf_l  ,ainarea*cf  _2) ;  */ 

printf ("\n\n\nMax  X-chord=  y,f  Max  Y-chord=  Y.f  Max  Area=  'it  ", 
xaax*cf _1 ,yaax*cf_l ,aaxarea*cf _2) ; 

printf ("\nMin  X-chord=  'it  Min  Y-chord=  y.f  Min  Area=  y.f  ", 
xain*cf _1 ,yain*cf_l ,ainarea*cf _2) ; 

/*  Close  open  file  and  free  allocated  aeaory  */ 

f  closed eatdata) ; 
freeCxptr) ; 
free(yptr) ; 
free(aptr) ; 
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APPENDIX  H:  PROGRAM  FILE:  speckle.c 


/♦PURPOSE  :  Provides  f or  calculating  the  speckle  index  of  an  image. 

The  resulting  value  is  used  as  a  Measure  to  evaluate  the 
sucess  of  filtering  an  image  to  reduce  the  speckle  noise. 
This  routine  is  called  in  all  filter  algorithms. 

— >  OPERATOR  must  enter  a  lumber  of  ROWS  and  COLUMIS 
THE  PROGRAM  uses  to  calclate  the  speckle  index. 

ROWS  - >  Maximum  =  480  Minimum  =  1 

COLS  - >  Maximum  =  512  Minimum  =  1 

♦/ 

•include  "thesis. h" 
speckle () 

{ 

extern  int  srow,  scol; 

int  m,  n.  x,  y,  nn; 

long  smax,  smin; 

long  sdata[NUM] ; 

unsigned  long  deviation,  ssum; 

float  smean,  slocal,  stotal,  tot; 
float  spkl index; 

if(  srow  >=  HROW  )  srow  =  IROW; 

if (  scol  >=  NCOL  )  scol  =  ICOL; 

tot  =  (long)srow  ♦  scol; 

/♦  Commence  calculation  ♦/ 

printf ("\n\n\tCalculating  Speckle  IIDEX ..."); 

for  (  n  =  1  ;  n  <  srov-2  ;  n++) 

{ 

for  (  m  =  1  ;  m  <  scol-2  ;  m+  +  ) 

{ 

smax  =  0; 

smin  =  260; 

ssum  =  0; 

nn  =  0; 


(ii 


1 or(  y  =  (  n-1  )  ;  y  <  (n+2)  ;  y++  ) 

{ 

for(  x  =  (  a-1  )  ;  x  <  (  a+2  )  ;  x++) 

{ 

■data [as]  =  rpixel(  x,y  ); 

■sum  +=  (long)sdataCnn] ; 

if(  saax  <  «data[nn]  )  taax  =  sdata[nn] ; 

if(  sain  >  •dataCnn]  )  sain  =  adataCnn]; 

nn++; 

> 

> 

deviation  =  saax  -  sain; 
il (  b sub  ==  0  )  saean  =  1; 

■a«an  =  ssun  /  IUH  ; 

■local  =  (float)deviation  /  saean; 

■total  +=  slocal; 

} 

> 

spklindex  =  stotal  /  tot; 

printf  ("\n\n\tThe  Calculated  '/.d  by  %d  speckle  index  is  */.f  ", srow,»col, spklindex ) ; 

> 

/* .  EID  Speckle  Function  . */ 
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APPENDIX  I:  PROGRAM  FILE:  vir.arry.c 


/♦PURPOSE  :  These  functions  provide  f or  ths  setup  of  a  disk  drive  virtual 
array  that  can  be  indexed  as  if  it  were  in  the  calling 
programs  data  storage  area.  The  routines  provide  for  disk 
file  access  when  required  to  retrieve  data  elements  fro*  the 
array.  This  allows  the  size  of  a  declared  array  to  be  limited 
only  by  the  amount  of  DISK  SPACE  available.  See  "thesis. h" 
header  file  for  setup  and  MACRO  routines  that  support  these 
functions.  Also  program  lstat.c  fully  implements  the  array 
routines . 

*/ 

Winclude  "thesis. h" 

/*  Virtual  Array  Access  Routines  ♦  / 

int  init_v_array (f ilename ,rec_size,f ilchar) 

char  ♦filenasie,  f ilchar; 

int  rec_size; 

{ 

long  size; 

FILE  ef; 

f  =  fopen(f ilename, "wb") ; 
if  (f  ! =  IULL) 

{ 

size  =  0; 

fwrite(Asize,4, 1 ,f ) ; 
fwrite(Arec_size,2, 1 ,f ) ; 
f write (tf ilchar , 1 , 1 , f ) ; 
f close(f ) ; 
return(l) ; 

> 

else 

return(IIULL) ; 

> 

VACB  *open_v_array (f ilenaae , buf f er_s ize) 

char  ♦filename; 

int  buffer.size; 

VACB  ev.array; 

char  ebuf.ptr; 

int  i ; 

char  f ilchar; 


/*  write  array  size  of  zero  ♦/ 
/*  and  array  element  size  ♦/ 
/♦  and  fill  char  ♦/ 
/♦  to  file  header  */ 


/*  allocate  virtual  array  control  block  */ 


v.array  =  (VACB  *)  malloc(sizeol (VACB) ) ; 
il  (v_arTay  ==  BULL)  re turn (BULL) ; 

/*  open  virtual  array  file  */ 

v_array->f ile  =  lopen(lilename,"r+b"); 
il  (v_array->l ile  ==  BULL) 

{ 

lree(v_array) ; 
return(BULL) ; 

>; 


/*  get  array  size  and  element  size  lor  control  block  */ 

lread(*v_array->size,4, 1 ,v_array->lile) ; 
lread(*v_array->elsize,2, 1 ,v_array->lile) ; 
lread(*lilchar , 1 ,1 , v_array->lile) ; 
v_array->bul_elsize  =  v_array->elsize  +  4; 

/*  allocate  buller  •/ 


v_array->buller  =  (char  *)  malloc(v_array->bul_elsize  *  (buller_size  +  1)) 
il  (v_array->bulfer  ==  BULL) 

{ 

lclose(v_array->lile) ; 
lree(v.array) ; 
return(BULL) ; 

>; 

v_array->bul_size  =  buller_size; 

/*  set  up  blank.rec  using  the  lill  character  in  array  header  */ 

/*  lor  initializing  new  array  elements  */ 

bul.ptr  =  v_array->buller  +  v_array->bul_elsize  *  v_array->bul_size; 

v_array->blank_rec  =  bul_ptr  +  4; 

lor  (i  ~  0;  i  <  v_array->bul_elsize;  i++) 

*bul_ptr++  =  lilchar; 

/*  set  record  index  negative  lor  all  buller  elements  */ 

buf_ptr  =  v_array->buller ; 

lor  (i  =  0;  i  <  v_array->bul_size ;  i++) 

< 

•((long  *)bul_ptr)  =  -IL; 
bul_ptr  +=  v  array->bul  elsize; 

>; 

return (v.array) ; 
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> 


void 

VACB 

{ 

int 

char 

long 


close.v.array(v.array) 

•v.array; 


i; 

•buf.ptr ; 

rec.index,  file.offset; 


buf.ptr  =  v_array->buffer; 


/*  Hush  buffer  */ 

for  (i  =  0;  i  <  v_array->buf_size;  i++) 

•C 

/*  check  each  element  index  */ 

/*  if  element  present;  write  it  to  disk  */ 

rec_index  =  *((long  *)buf_ptr); 
if  (rec_index  >=  0) 

{ 

file.offset  =  header  +  rec.index  *  v_array->elsize; 
f seek( v_array->f ile , f ile_of f set , 0) ; 

fwrite(buf_ptr  +  4,  v_array->elsize , 1 ,  v_array->f ile) ; 

>; 

buf_ptr  +=  v_array->buf  elsize; 

>; 

free(v_array->buff er) ;  /*  de-allocate  buffer  */ 

f close(v_array->f ile) ;  /*  close  array  file  */ 

free(v_array) ;  /*  de-allocate  VACB  */ 


void  *access_v_rec(v_array , index) 

VACB  *v_array; 

long  index; 

{ 


char  *buf_ptr; 

int  buf .index; 

long  rec.index,  temp.index; 


/*  calculate  buffer  address  of  referenced  element  */ 
buf.index  =  index  '/,  v_array->buf_size; 

buf.ptr  =  v_array->buf f er  +  buf.index  *  v_array->buf_elsize; 
rec.index  =  ‘(long  *)buf_ptr; 
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/*  if  element  present,  return  its  buffer  address  ♦/ 
if  (rec.index  ==  index)  return(buf _ptr  ♦  4) ; 

/♦  if  element  doesn’t  exist,  extend  the  file  ♦/ 

if  (index  >=  v_array->size) 

f seek (v_array->f ile, 0,2) ; 

for  (temp.index  =  v_array->size;  teep_index++  <=  index;  ) 

fHrite(v_array->blank_rec,  v_array->elsize,  1,  v_array->f ile) ; 
v_array->size  =  index  ♦  1; 
f seek (v_array->f ile. 0,0) ; 

fwrite(kv_array->size,  4,  1,  »_array->f ile) ; 

>; 

/*  if  buffer  slot  is  occupied  by  another  elenent,  ♦/ 

/♦  save  it  to  disk  ♦/ 

if  (rec.index  >=  0) 

{ 

fseek(v_array->f ile,  rec.index  *  v_array->elsize  +  header,  0); 
fwrite(buf_ptr  +  4,  v_array->elsize,  1,  v_array->f ile) ; 

>; 

/♦  read  referenced  element  into  buffer  slot  ♦/ 

f seek(v_array->f ile,  index  *  v_array->elsize  +  header,  0); 
fread(buf_ptr  +  4,  v_array->elsize,  1,  v_array->f ile) ; 

♦((long  *)buf _ptr)  =  index; 

/♦  return  address  of  element  in  buffer  ♦/ 


} 

/*... 


*/ 


return(buf _ptr  +  4); 


APPENDIX  J:  PROGRAM  FILE:  genfunc.c 


/♦PURPOSE  :  Those  general  functions  provide  f or  routine  evolutions 
that  occur  a  number  of  times  in  the  Image  processing 
routines  developed  for  analysis  of  speckle  reduction 
algorithms.  Some  require  support  of  other  functions 
and  are  not  totally  independent. 

♦/ 

#include  "thesis.h" 

/*  initial  ITEXPC  Board  setup  */ 

start it () 

{ 

sethdw ( REGB ASE , MEKB ASE , MEMORY ) ; 
setdim(XSIZE,YSIZE, DEPTH) ; 
aclear ( IXS , I YS , ICOL , IROW , COLOR) ; 
fgon() ; 
initial ize() ; 
clsO ; 

select_mem(MEH_A) ; 
display _mem(MEM_A) ; 

> 

/*  Get  integer  keyboard  input  */ 

get i( name, iptr) 
char  ♦name ; 
int  ♦iptr; 

{ 

printf("  */.s  ".name); 
scanf("  y.d",iptr); 

} 


/*  Get  floating  point  keyboard  input  ♦/ 

getf (rname , riptr ) 
char  ♦rname; 
int  eriptr; 

{ 

printf('"/,s  ".rname); 
scanf("  */,f"  .riptr) ; 

} 


/*  Pause  routine  used  to  prompt  for  operator  intervention*/ 
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uait() 

{ 

print! ("\n\n  Press  Return  to  continue"); 
!llush(stdin) ; 
get char () ; 

> 


/*  returns  minimum  integer  value  from  input  array  passed  •/ 

getmin (data, nun) 

int  *data; 
int  nun; 

int  t,  min; 

for  (nin=  data[0],  t=l  ;  t<nun  ;  t++) 
il(  dataCt]  <  nin  )  nin  =  data[t] ; 

return  nin; 

> 

/*  returns  naxinun  integer  value  iron  input  array  passed  */ 

getnax (data, nun) 

int  *data; 
int  nun; 

{ 

int  t,  nax; 

for  (nax=  data[0] ,  t=l  ;  t<nun  ;  t++) 
if(  data[t]  >  nax  )  nax  =  dataCt]; 

return  nax; 

> 

/*  DOS  Systen  Call  to  clear  Display  Screen  */ 

cls() 

{ 

_clearscreen(_GCLEARSCREEH) ; 

> 

/*  routine  to  read  an  inage  fron  disk  */ 
readitQ 

/♦progran  uses  the  following  externals:  filename, coniine*/ 

{ 

extern  char  filename [20] ,  coniine [200] ; 
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int 

int 


errval ; 
ch; 


whiled) 

{ 

printl ( " \n\nEITER  IMAGE  FILEIAME  as  (DEV : name. IMG) ->") ; 

ff lush(stdin) ; 

scanf  ("'/.s"  .filename) ; 

errval  =  readim(IXS, IYS.ICOL.IROW, filename, comline) ; 
if (errval  ==  0)  { 

printl ("\nFILE->  ’/.s  \n\nC0MMEIT :  Xs" .filename, comline) ; 
break; 

} 

if( errval  !=  0)  { 

switch(errval) 

{ 

case  FILE.ERROR: 

print! ("Error  opening  file\n"); 
break; 

case  FORMAT. ERROR: 

print! ("Unknown  file  formatW); 
break ; 

case  READ.ERROR: 

printl("Error  Reading  lile\n"); 
break; 
default : 

print! ("?  Unknown  Error  ?Xd\n" .errval) ; 
break; 

> 

printl ("Try  Again??  ...Yes(y)  /  Io(n)\n"); 
ch  =  getch() ; 

if(ch  ==  *y»  ||  ch  ==  ’y’)  continue; 
if(ch  ==  *N’  II  ch  ==  'n’)  exit(0); 

> 

> 

> 


/*  routine  to  save  an  image  to  disk  */ 
saveitO 

/•program  uses  the  following  externals:  *f ilename , *comline*/ 

{ 

extern  char  f ilename[20] ,  comline [200] ; 

int  format,  errval; 

int  ch; 

initialized ; 
while(l) 

{ 


print! ("\n\nEITER  FILENAME. . . (DEV:name. IKG)->") ; 

f f lush(stdin) ; 

scan!  ("’/,s"  .filename) ; 

geti("\n  Enter  File  Compression  Value  ?. . . (0/1)" .Aformat) ; 
printf("\n\n  Enter  Image  Comment  ..to  continue..  (Max  200  CHAR)..."); 
fflush(stdin) ; 
gets (coaline) ; 

if (format  ==  l)printf ("\n\nStoring  Image  Using  Compression  —  Please  wait  M") 
if (format  ==  0)printf ("\nStoring  Image  —  Please  wait!!"); 

errval  =  saveim(IXS, IYS,IC0L,IR0V, format .filename, coaline) ; 

if (errval  ==  0)  { 

print! ("\nlmage  save  completed  Sat”); 
break; 

> 

if (errval  !=  0)  { 

print! ("\n\a  Error  saving  file!!"); 
if (errval  ==  ALLOCATIOI.ERROR)  { 

print! ("\nlnsuff icient  Disk  Space"); 

> 

if (errval  ==  VRITE_ERR0R)  { 

print! ("\nError  writing  file  or  values"); 

> 

printf ("\n\tTry  Again??  ...Yes(y)  /  Ho(n)"); 
ch  =  getchO  ; 

if(ch  ==  »Y*  ||  ch  ==  ’y’)  continue; 
if ( ch  ==  >»’  ||  ch  ==  ’n’)  break; 

> 

> 

> 

/*  counter  step  routine  */ 

step(u,v) 
int  *u,*v; 

{ 

♦u  +=  1; 

if (*u  >=  HIGH)  { 

*u  =  1; 

*v  +-  1; 

> 

> 

/*  threshold  routine  for  feat.id.c*/ 

thresh( ) 

{ 

int  ans ,c, option; 
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unsigned  int  val; 


while(l) 

< 

initlutsO ; 
initial ize() ; 

get i ( *' \n\n\t EITER  IEW  THRESHOLD  LEVEL  (0-255)  :M,*val); 
setlut(IIPUT,0) ; 
threshold(IIPUT,0, HIGHEST, val) ; 
setlut(RED.O) ; 

threshold(RED,0, HIGHEST. val); 
setlut (GREEI.O) ; 
threshold (GREEI , 0 , HIGHEST , val) ; 
setlut (BLUE, 0) ; 

threshold(BLUE , 0 .HIGHEST , val) ; 
setlut (IHPUT.O); 

printf  ( "\n\nSATISFIED  WITH  THRESHOLD  AT  ’/.d?. .  .Yes(y)  /  Io(n)  .  .  . "  ,  val) ; 
ans  =  getchO ; 

if (ans  ==  ’ Y ’  II  ans  ==  ’y 1 )  break; 

> 

■aplut (IIPUT , 0 , IXS , I YS , HCOL , IROV ) ; 
initialized ; 


/*  threshold  routine  for  threshit.c  */ 

threshitO 

{ 

int  ans ,c, opt ion; 

int  *val; 

while(l) 

{ 

initlutsO ; 
initialized ; 

ge t i ( " \n\nEITER  HEW  THRESHOLD  LEVEL  (0-255)  .".Aval); 

printf("%d\nM,val); 

setlut ( IIPUT .0) ; 

threshold( IIPUT, 0, HIGHEST, val) ; 

setlut(RED.O) ; 

threshold( RED , 0 , HIGHEST , val ) ; 
setlut (GREEK, 0) ; 
threshold(GREEK,0, HIGHEST, val) ; 
setlut (BLUE ,0) ; 

threshold(BLUE , 0 .HIGHEST , val ) ; 
setlut ( IIPUT, 0) ; 

printf  ("\n\nSATISFIED  WITH  THRESHOLD  AT  '/.d?.  .  .Yes(y)  /  Io(n)  .  . . " , val) ; 
ans  =  getchO ; 
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if (ans  ==  *Y»  | i  ans  ==  ’y')  break; 


> 

print! (”\n\t\t  0:  SAVE  THE  THRESHOLD  IMAGE?..."); 

printl("\n\t\t  1:  QUIT  (Leave  thresholded  image  !or  further  processing)...") 
print! ("\n\t\t  2:  RESTORE  system  to  original  input  image..."); 
geti(“\n\t  Select  option  by  HUMBER:  “.koption); 
avitch(option) 

case  0: 

maplut ( IIPUT , 0 , IXS , ITS , ICOL , IROW ) ; 
saveitO ; 
break; 
case  1: 

maplut ( IIPUT ,0,11s, IYS , ICOL , IROW ) ; 
break; 
case  2: 

initializeQ ;  /*initialize  without  mapping  operation*/ 

break; 
default : 
break; 

> 

> 

/*  function  to  compute  standard  deviation  of  image  */ 

dev(stddev) 

float  *stddev; 

float  var,  sum,  sqsum,  sqvalue,  dim,  value; 

long  X ,  y ; 

sum  =  sqsum  =  0; 

dim  =  (float) IROW  *  ICOL; 

for(  y  =  0  ;  y  <=  IROW-1  ;  y++) 

{ 

for(  x  =  0  ;  x  <=  ICOL-i  ;  x++) 

i 

value  =  rpixel(x.y); 
sqvalue  =  value  *  value; 
sum  +=  value/dim; 
sqsum  ♦=  sqvalue/dim; 

> 

> 

var  =  sqsum  -  (sum*sum); 

♦stddev  =  sqrt(var)  ; 


/*  function  to  compute  scale  factor  for  sizing  routine  */ 
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scale(f actor) 
float 


♦factor; 


int  inch,  pixel,  flag; 

float  dinch,  dpixel;  /*  default  length  values  ♦  / 

/♦  dinch  and  dpixel  are  to  be  Modified  accordingly  ♦/ 
dinch  =  1.00; 

dpixel  =  BO. 00; 

printf ("\n\tUtilizing  standard  SCALE_FACT0R? . . . Yes  (y)  /  lo  (n)"); 
flag  =  getchO; 

if (flag  ==  »Y*  I i  flag  ==  *y') 

{ 

♦factor  =  ( (f loat)SCALE_FACTOR) ; 

> 

if (flag  ==  ’I’  I  I  flag  ==  'n') 

{ 

printf ("\n\tPress  ’0’  if  using  default  values  for  length  in  inches  or  pixels"); 
geti("\n\tEnter  the  length  of  object  in  inches (default) :",Ainch) ; 
geti("\n\tEnter  the  length  of  object  in  pixels (def ault) :",*pixel) ; 

if (inch  ==  0) 

{ 

inch  =  dinch; 

> 

if(pixel  ==  0) 

{ 

pixel  =  dpixel; 

> 

♦factor  =  ( (float)inch)  /  ( (float)pixel) ; 

> 

printf ("\n\tSCALE_FACT0R  is  */,f \n" , ♦factor) ; 

> 
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APPENDIX  K:  C  COMPILER  OPERATIONS 

The  following  batch  files  list  the  C  compiler  directives  used  to  compile  and  link 
all  C  programs.  The  C  compiler  files  ACCOMP.BAT  and  DACCOMP.BAT  are 
listed  below.  ACCOMP.BAT  contains  the  switches  necessary  to  utilize  CODEVIEW 
debugger,  and  ACCOMP2.BAT  is  for  final  program  optimization.  The  F  drive  listed 
in  the  PATH  statement  of  these  batch  files  is  due  to  the  fact  that  the  programs  were 
stored  on  a  Bernoulli  disk. 


ACCOMP.BAT: 

path  c : \ ; c: \dos; c:\bs jet ;f :\k«dit.*;c:\Bousel;c:\xtree; c:\lsrctrl; . . . 
f :\;i :\«sc\bin;f :\msc\tmp;f :\msc\src;l :\*sc\wrk;f :\a8c\lib; . . . 
f:\nsc\include 
set  include=\B8c\include 
set  lib=\msc\lib;\pcplus\it«x\lib 
set  tmp=\msc\tmp 

cl  /c  /Fs  /Zi  /GO  /AL  /Fm  */.l.c 

ACC0MP2.BAT 

path  c:\;c:\dos;c:\*isjet;f  :\kedit c:\mousel ;c : \xtree; c: MsrctzT ; .  .  . 
f :\;f :\®sc\bin;f :\msc\tmp;f :\msc\src;f :\msc\wrk;T :\msc\lib; . . . 
T:\msc\include 
set  include=\msc\include 
set  lib=\msc\lib ; \pcplus\itex\lib 
set  tmp=\msc\tmp 

cl  /c  /Fs  /Ze  /GO  /AL  */.! .  c 


The  following  batch  files,  LTHESIS.BAT  and  DTHESIS.BAT,  link  the  programs 
to  the  libraries  and  supporting  function  files.  The  speckle. c  should  not  be  included 
within  these  files  when  linking  the  programs  sizeit.c  and  feat-id. c  due  to  an  undefined 
external  problem  associated  with  the  external  variables  SROW  and  SCOL.  LTHE¬ 
SIS.BAT  is  used  for  debugging  applications,  and  DTHESIS.BAT  is  for  final  program 
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optimization.  It  is  also  unnecessary  to  link  feat-id. c  and  sizeit.c  to  the  virtual  array 
functions  contained  within  vir.arry.c. 


LTHESIS.BAT: 

IF  *0T  ERRORLEVEL  1  LIIK  /IOD  /CO  /E  Xl  genfunc  speckle  vir.arry  . . . 
, , , ITEXPCML . LIB  LLIBCE . LIB ; 


DTHESIS.BAT: 

IF  *0T  ERRORLEVEL  1  LIIK  /IOD  ,?  */.l  genfunc  speckle  vir.arry... 

, , .ITEXPCML. LIB  LLIBCE. LIB; 

To  create  the  menu  format  for  program  display,  a  C  program  named  hologram. c 
was  generated.  To  link  all  image  processing  programs  to  this  module,  all  program 
umain()r’  statements  were  modified  so  that  hologram. c  was  made  the  main  program. 
All  image  processing  programs  were  then  linked  to  hologram. c  after  compilation.  The 
ACCOMP2.BAT  file  was  used  for  compilation,  and  the  PTHESIS.BAT  file  was  used 
for  subsequent  module  linking.  Due  to  the  length  of  the  LINK  statement,  the  user 
must  enter  the  libraries  ITEXPCML  and  LLIBCE  manually  into  the  keyboard  each 
time  the  LINK  process  is  accomplished. 

PTHESIS.BAT: 

IF  NOT  ERRORLEVEL  1  LINK  /IOD  /E  ’/.I  table  speckle  vir.arry... 

genlunc  threshit  sizeit  feat.id  geofil. . 
lstat  2sigma 
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